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ABSTRACT 


This  report  presents  design  procedures  and  parametric  that  can  be  used 

to  avoid  flutter  and  divergence  of  hydrofoil  strut-foil  systems.  Flutter  stability  of 
conventional  T-foils  can  be  achieved  by  scaling  existing  prototype  systems  which 
are  stable;  scaled  flutter  models  must  be  used  to  determine  the  stability  of  n  w 
designs.  Several  stable  T-foil  designs  are  available  for  craft  with  displacements  of 
up  to  300  tons  at  subcavitating  speeds  (below  50  knots).  The  availability  of 
higher  speed  designs  for  ships  of  practical  size  is  limited  to  one  system  that  is 
stable  to  62  knots.  Divergence  stability  of  subcavitating  systems  can  be  ensured 
by  theoretical  analysis  using  an  available  computer  program. 


ADMINISTRATIVE  INFORMATION 

This  work  was  authorized  and  funded  under  the  Hydrofoil  Development  Program  of  the 
Naval  Ship  Systems  Command,  Subproject  S4606,  Task  1703,  Work  Unit  1 153-003. 
Additional  support  was  given  by  the  Naval  Material  Command  Program  Element  62754N, 
Task  Area  ZF43421001,  Work  Unit  1520-001. 

INTRODUCTION 

In  this  report,  results  of  flutter  and  divergence  research  have  been  rephrased  for 
practical  application  by  the  designer.  Procedures  are  included  both  for  producing  a  pre¬ 
liminary  design  and  for  evaluating  a  final  design. 

The  appendages  of  high-speed  ships  must  be  free  of  hydroelastic  instability  (flutter  and 
divergence).  To  date,  all  such  appendages  have  been  free  of  these  problems.  Divergence 
stability  was  achieved  by  using  calculations  that  were  sufficiently  accurate  for  design 
purposes.  However,  because  flutter  boundaries  could  not  be  calculated  accurately,  other 
design  criteria  were  responsible  for  the  flutter  stability  of  existing  struts  and  foils.  Those 
were  based  on  maximum  expected  stress  leveis.  Specifically,  struts  and  foils  suitable  foi  use 
in  rough  seas  were  rigid  enough  to  resist  flutter  and  divergence.  Despite  the  proven  stability 
of  existing  designs,  hydroelastic  prediction  techniques  are  needed  for  reassurance  on  the 
stability  of  future  conventional  appendage  designs  as  well  as  for  the  successful  design  of 
unconventional  appendages. 

Notwithstanding  recent  significant  advances  in  understanding  flutter  and  divergence, 
neither  existing  data  nor  theoretical  analyses  are  presently  adequate  to  offer  accurate  pre¬ 
dictions  of  flutter  boundaries.  Therefore,  the  designer  of  future  systems  is  offered  only  an 
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This  report  presents  design  procedures  and  parametric  trends  that  can  be  used 
to  avoid  flutter  and  divergence  of  hydrofoil  strut-foil  systems.  Flutter  stability  of 
conventional  T-foils  can  be  achieved  by  scaling  existing  prototype  systems  whkh 
are  stable;  scaled  flutter  models  must  be  used  to  determine  the  stability  of  new 
designs.  Several  stable  T-foil  designs  are  available  for  craft  with  displacements  of 
up  to  300  tons  at  subcavitating  speeds  (below  50  knots).  The  availability  of 
higher  speed  designs  for  ships  of  practical  size  is  limited  to  one  sy^em  that  is 
stable  to  62  knots.  Divergence  stability  of  subcavitating  systems  can  be  ensuied 
by  theoretical  analysis  using  an  available  computer  program. 
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INTRODUCTIOi  1 

In  this  report,  results  of  flutter  and  divergence  research  have  been  rephrased  for 
practical  application  by  the  designer.  Procedures  are  included  both  for  producing  .t  pre¬ 
liminary  design  and  for  evaluating  a  final  design. 

The  appendages  of  high-speed  ships  must  be  free  of  hydroelastic  instability  (flutter  and 
divergence).  To  date,  all  such  appendages  have  been  free  of  these  problems.  Divergence 
stability  was  achieved  by  using  calculations  that  were  sufficiently  accurate  for  design 
purposes.  However,  because  flutter  boundaries  could  not  be  calculated  accurately,  other 
design  criteria  were  responsible  for  the  flutter  stability  of  existing  struts  and  foils.  These 
were  based  on  maximum  expected  stress  levels.  Specifically,  struts  and  foils  suitable  tor  use 
in  rough  seas  were  rigid  enough  to  resist  flutter  and  divergence.  Despite  the  proven  stability 
of  existing  designs,  hydroelastic  prediction  techniques  are  needed  for  reassurance  on  the 
stability  of  future  conventional  appendage  designs  as  well  as  for  the  successful  design  of 
unconventional  appendages. 

Notwithstanding  recent  significant  advances  in  understanding  flutter  and  divergence, 
neither  existing  data  nor  theoretical  analyses  are  presently  adequate  to  offer  accurate  pre¬ 
dictions  of  flutter  boundaries.  Therefore,  the  designer  of  future  systems  is  offered  only  an 
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indirect  approach,  utilizing  either  the  performance  of  existing  full-scale  hydrofoil  structures 
or  simulation  with  scaled  flutter  models.  Quantitative  stability  boundaries  can  be  established 
on  the  basis  of  the  general  scaling  laws  provided. 

The  present  report  gives  a  method  of  theoretical  analysis  which  is  sufficiently  accurate 
for  divergence  prediction  and,  in  some  cases,  flutter  prediction  also.  Subcavitating  and  cavi- 
tating  systems  are  treated  separately.  The  report  includes  a  brief  description  of  flap  reversal, 
a  hydroelastic  effect  which  nullifies  the  control  effect  of  a  flap  and  therefore  must  be 
avoided.  It  is  not  an  instability,  however.  A  typical  inverted-T  strut-pod-foil  system  is 
shown  in  Figure  1.  Most  available  research,  however,  deals  with  struts  without  foils.  On 
hydrofoil  craft,  T-foils  appear  to  be  relatively  more  vulnerable  to  hydroelastic  instability  than 
do  inverted-pi  foils,  the  other  common  configuration.  Inverted-pi  foils  are  more  stable  be¬ 
cause  the  struts  and  inboard  foil  sections  are  multiply  constrained. 


DESCRIPTION  OF  FLUTTER,  DIVERGENCE,  AND  FLAP  REVERSAL 
FLUTTER 

Flutter  is  the  unstable  oscillation  of  a  structure  composed  of  lifting  surfaces  traveling 
through  a  fluid.  The  oscillation  increases  exponentially  in  amplitude  until  the  structure 
yields  or  the  flow  changes  in  some  way.  The  strut-foil  systems  on  hydrofoil  craft  and  the 
strut-like  rudders  expected  to  be  used  on  surface  effect  ships  are  examples  of  structures  which 
may  be  vulnerable  to  flutter. 

Flutter  characteristics  of  structures  are  specified  in  terms  of  speed  boundaries  which 
separate  regions  of  posirive  and  negative  damping.  The  boundaries  are  composed  of  critical 
flutter  speeds,  also  known  as  flutter  inception  speeds,  and  correspond  to  zero  damping  of 
the  system.  At  a  critical  flutter  speed,  an  oscillation  triggered  by  a  deflection  of  the 
structure  will  continue  at  constant  amplitude.  Above  this  speed,  an  oscillation  will  increase 
in  amplitude.  It  is  expected  that  there  will  be  sufficient  unsteadiness  :.n  water  flow  to  trigger 
oscillation  at  or  above  the  flutter  speed. 


Bending  and  Torsional  Flutter  Modes 

Experimental  and  theoretical  results  indicate  that  cantilevered  struts  undergo  flutter  in 
two  different  hydroelastic  modes.  Based  on  the  predominant  structural  mode  shapes 
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involved  these  modes  have  been  designated  as  bending  flutter  and  torsional  flutter.1  Since 
the  flutter  characteristics  of  each  mode  are  significantly  different,  it  is  important  to  identify 
the  mode  to  be  expected  for  a  given  strut. 

The  flutter  mode  of  a  strut  can  be  determined  by  examining  the  mode  shape  of  its 
second  vibration  mode  in  water,  except  in  a  transition  region  where  strong  coupling  of 
structural  modes  occurs.  Mode  shapes  are  designated  by  their  similarity  to  the  uncoupled 
mode  shapes  of  a  cantilever  beam,  which  are  shown  in  rigure  2.  All  struts  exhibit  a  funda¬ 
mental  (lowest  frequency)  vibration  mode  shape  tha'  resembles  first  bending.  However,  the 
marked  difference  exhibited  by  struts  in  their  second  (next-to-lowest  frequency)  mode 
shapes  permits  identification  of  the  flutter  mode.  If  the  second  mode  shape  is  predominantly 
second  bending,  the  strut  will  undergo  flutter  in  a  predominantly  first  bending  mode  shape. 
This  type  of  flutter  is  called  bending  flutter,  and  the  struts  are  known  as  bending-tyne 
struts.  If  the  second  mode  shape  of  the  strut  is  predominantly  first  torsion,  the  strut  v'M 
undergo  flutter  in  a  predominantly  first  torsion  mode  shape.  This  type  of  flutter  is  called 
torsional  flutter,  and  the  struts  are  known  as  torsion-type  struts. 

Strut*  with  little  or  no  tip  weighting  are  usually  bending-type  struts,  and  those  with 
relatively  heavy  pods  are  usually  torsion-type  struts.  There  is  a  transition  region  in  which 
the  second  vibration  mode  of  a  strut  has  strong  components  of  both  second  bending  and 
first  torsion,  with  neither  predominating.  Struts  in  this  region  have  moderately  weighted 
pods  or  medium  to  large  foiis.  It  is  impossible  to  deduce  which  flutter  mode  will  occur  in 
this  transition  region. 

Typical  hydroelastic  mode  characteristics  for  bending-type  and  torsion-type  struts  are 
shown  in  Figures  3  and  4,  respectively.  These  drawings  were  adapted  from  published  calcu¬ 
lations  and  experimental  results.1 

As  shown  in  Figure  3,  bending  flutter  corresponds  to  an  instability  in  the  “new  mode,” 
a  mode  which  does  not  exist  at  low  speed.  The  new  mode  first  occurs,  highly  damped,  at 
intermediate  speed  and  decreases  rapidly  in  damping  just  below  the  critical  flutter  speed  Ur. 
The  new  mode  has  a  first  bending  mode  shape.  Its  frequency  remains  very  low  and  close  to 
that  of  Mode  1.  Because  of  its  damping  behavior,  the  new  mode  cannot  be  detected  in 
practice  until  shortly  before  flutter  is  encountered. 

An  entirely  different  hydroelastic  mode  is  responsible  for  torsional  flutter,  as  shown  in 
Figure  4.  The  damping  of  Mode  2  remains  low  throughout  the  subcritical  speed  range, 
rising  on.  •  slightly  above  its  zero  speed  value  before  becoming  negative.  Flutter  occurs  in  a 
predominantly  first  torsion  mode  shape. 


Besch,  P.K.  and  Y.-N.  Liu,  “Bending  Flutter  and  Torsional  Flutter  of  Flexible  Hydrofoil  Struts,"  Ninth  Symposium  on 
Naval  Hydrodynamics,  Paris,  France  (Aug  1972);  also  NSRDC  Report  4012  (Feb  1971).  A  complete  listing  of  references  is 
given  on  pages  73-76. 


4 


Figure  2  -  Uncoupled  Mode  Shapes  of  Cantilev  :r  Beam 


SPEED  U 


SPEED  U 

Figure  3  -  Typical  Hydroelastic  Mode  Characteristics  for  Bending- 

Type  Struts 


FREQUENCY 


Figure  4  -  Typical  Hydroclastic  Mode  Characteristics  for  Torsion- 

Type  Struts 


Calculations  imply  that  at  very  high  speeds  a  "new  mode”  appears  in  the  torsional 
flutter  region.  As  strut  characteristics  vary  toward  the  bending  region,  the  stability  of 
Mode  2  increases  and  that  of  the  new  mode  decreases.  The  transition  from  torsional  flutter 
to  bending  flutter  appears  to  involve  a  crossover  of  the  critical  flutter  boundaries  of  these  two 
modes. 

Strut  system  characteristics  are  composed  of  inertial,  elastic,  and  damping  characteristics. 
The  change  from  bending  flutter  to  torsional  flutter  occurs  as  a  result  of  inertial  and  elastic 
changes  in  the  strut  system.  Damping  is  not  usually  considered  to  be  variable.  Inertial 
changes  have  a  complicated  effect  on  flutter  speed;  these  are  qualitatively  well  understood 
for  strut-pod  systems  and  less  so  for  T-foiis  Effects  of  elastic  changes  are  only  approximately 
known. 


Effect  of  Changing  System  Inertia 


To  some  extent,  important  system  inertial  properties  can  be  summarized  in  the  parameter 
generalized  mass  ratio  p.  This  parameter  is  the  rafio  of  structural  and  fluid  inertia  for  a 
given  motion  of  the  strut-fluid  system.  Generalized  mass  ratios  have  been  calculated  for 
first  bending  modes  of  a  T-foil  by  Mitchell  and  Rauch.2  To  permit  a  simple  calculation, 
approximations  have  been  derived  by  assuming  pure  bending  motion  and  pure  torsional 
motion  to  simulate  bending  flutter  and  torsional  flutter,  respectively.1  For  bending  motion 
of  strut-pod  systems, 


Mstrut  +  Mpod 

^bending  “  M« -  +  M* 

strut  pod 


where  M*(Ut  =  irpc2  C/4. 

For  torsional  motion  of  strut-pod  systems, 


^y ,  strut  *  ^y,  pod 

^torsion  "  r#  +  I# 

y’  strut  y,  pod 


*y,  strut 


jr/>c4r  I  I  /  ^xea  \Z  \ 

*  -wT  It  *(t 


2 

Mitchell,  L.  and  F.J.  Rauch,  Jr.,  "Dynamic  Tests  of  the  1/4-Scale  Models  of  the  80-Knot  Transiting  Strut-Foil  Systems 
for  the  FRESH  1  Hydrofoil  Test  Craft,”  Grumman  Aircraft  Engineering  Corporation,  Contract  NObs  86826  (Aug  1964). 
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The  expressions  given  for  strut  added  mass  M*rut  and  added  moment  of  inertia  I*  mu|  apply 
to  struts  of  uniform  chord  c  and  elastic  axis  location  *  Pod  added  mass  M*od  and  added 
moment  of  inertia  If  pod  were  approximated  in  the  present  work  by  values  given  by 
Kennard3  for  prolate  spheroids  having  lengths  and  maximum  diameters  equal  to  those  of  the 
pods. 

Model  results  indicate  that  bending-type  struts  with  solid  cross  sec.ions  undergo  bending 
flutter  for  values  of  Mbendbu  between  0.1  and  0.66.  Flutter  speeds  for  this  mode  of  flutter 
increase  monotonically  as  jubendint!  increases.  Flutter  frequencies  are  low;  reduced  frequencies 
knf  are  n<  ir  0.1.  Torsion-type  struts  with  solid  sections  undergo  torsional  flutter  for  values 
of  /i,orsion  between  0.66  and  6.2.  Flutter  speeds  decrease  rapidly  as  A*torSion  increases  above 
0.66,  and  they  may  reach  a  minimum  at  some  higher  value  of  MlorJion-  Reduced  frequencies 
knf  are  near  1.0.  Flutter  speeds  are  generally  lower  and  flutter  frequencies  higher  than  for 
bending-type  struts.  These  flutter  speed  relationships  are  shown  in  Figure  5. 

By  implication  from  Figure  5,  the  most  stable  inertial  configuration  of  a  strut  system  is 
that  which  places  the  strut  on  the  boundary  between  bending  flutter  and  torsional  flutter. 
This  condition  occurs  when  the  first  torsion  and  second  bending  frequencies  are  equal.  Con¬ 
sequently,  the  flutter  speed  of  a  strut  with  a  relatively  heavy  pod  can  be  raised  by  lightening 
the  pod  or  by  attaching  a  high  aspect  ratio  foil. 


Effect  of  Changing  System  Stiffness 

System  elastic  properties  determine  the  value  of  mass  ratio  at  which  the  bending  and 
torsional  flutter  boundaries  intersect  as  well  as  the  magnitude  of  flutter  speed  at  a  given  n. 

For  thin,  solid-section  struts  with  blunt-based  profiles,  the  intersection  between  bending  and 
torsional  flutter  occurs  at  n  =  0.66.  Solid  struts  with  airfoil-shaped  profiles  or  struts  con¬ 
structed  of  spars,  ribs,  and  skin  are  expected  to  have  approximately  the  same  flutter  boundary 
intersection,  but  these  have  not  been  investigated  as  thoroughly.  This  intersection  appears 
to  be  affected  by  elastic  axis  location  and  the  ratio  of  El  to  GJ.  For  otherwise  similar 
systems,  flutter  speed  magnitudes  are  related  to  the  magnitudes  of  El  and  GJ. 


3 

Kennard,  E.H.,  “Irrotational  Flow  of  Frictionless  Fluids,  Mostly  of  Invariable  Density,”  David  Taylor  Model  Basin 
Report  2299,  pp.  390-392  (Feb  1967). 


9 


BENDING 

FLUTTER 

REGION 


Figure  5  -  Variation  of  Flutter  Speed  Uf  as  a  Function  of  Generalized 
Mass  Ratio  ju  for  Simple  Struts  and  Struts  with  Pods 

(Value  of  miu  ratio  at  intersection  of  bending  and  torsional  flutter  modes  is  determined  by 
elastic  and  inertial  characteristics  of  structure;  0.66  occurs  for  solid  section  struts  with  thin, 
blunt-based  profiles.) 


Effect  )f  Changing  Structural  Damping 


Structural  damping  gives  a  vertical  offset  to  the  hydrodynamic  damping  curve,  thereby 
affecting  the  resulting  flutter  inception  speed  at  which  total  damping  becomes  zero.  The 
unstable  hydroelastic  mode  in  bending  flutter  has  a  rapid  decrease  in  damping  near  flutter 
and,  therefore,  is  little  affected  by  changes  in  structural  damping.  In  contrast,  torsional 
flutter  speeds  are  sometimes  very  sensitive  to  changes  in  structural  damping  because  hydro- 
dynamic  damping  is  low  at  all  speeds.  The  effect  of  adding  a  measured  value  of  structural 
damping  to  the  hydrodynamic  damping  predicted  for  the  torsional  flutter  mode  is  shown  in 
Figure  4. 

Flutter  has  been  obtained  experimentally  for  a  large  number  of  strut  models  at  speeds 
ranging  from  5.9  to  1 16  knots.  Most  of  the  data  are  for  struts  with  blunt-based  profiles. 
These  struts  appear  to  undergo  flutter  in  the  same  mode  as  struts  with  subcavitating  airfoil¬ 
shaped  profiles  but  at  different  speeds. 


DIVERGENCE 

Divergence  is  static  instability  of  a  lifting-surface  structure  traveling  through  a  fluid.  It 
occurs  when  steady  elastic  deformations  of  a  structure  lead  to  increases  in  hydrodynamic 

loading.  The  divergence  speed  corresponds  to  the  condition  in  which  the  increase  in  hydro- 

» 

dynamic  loading  is  exactly  balanced  by  the  increase  in  elastic  force.  Theoretically,  infinite 
deflections  occur  at  higher  speeds.  In  practice,  structures  yield  well  before  the  divergence 
speed  is  reached. 

There  are  four  common  misconceptions  about  divergence.  First,  divergence  is  not 
simply  a  “strength  problem.”  Any  failure  due  to  appicaching  divergence  can  be  avoided  by 
sufficiently  increasing  the  stiffness  of  the  structure,  but  the  failure  is  a  result  of  elastic 
deflections  which  change  the  loading  on  the  structure.  Any  design  procedure  must  include 
the  elastic  response  of  the  structure  to  the  applied  loading.  This  can  be  done  either 
iter? lively  or,  more  efficiently,  by  making  a  stability  calculation  similar  to  a  buckling  analysis. 
The  stability  boundary  is  independent  of  initial  angle  of  attack  although  angle  of  attack  will 
affect  the  speed  at  which  yield  occurs. 

Second,  structures  cannot  be  operated  at  speeds  ‘  up  to”  the  divergence  speed.  When  a 
divergence  instability  is  present,  deflections  increase  extremely  rapidly  over  the  last  one- 
third  of  the  speed  range  below  the  divergence  speed.  The  large  deflections  will  hamper 
operational  characteristics  of  hydrofoil  struts,  will  decrease  fatigue  lifetime  when  oscillatory 
loading  is  involved,  and  will  cause  yielding  well  below  the  divergence  speed. 
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Third,  struts  cannot  be  assumed  to  have  an  aerodynamic  center  location  at  quarter- 
chord.  This  assumption  has  led  designers  to  conclude  that  the  divergence  speed  of  a  strut  was 
infinite  when  actually  it  was  finite.  Divergence  analysis  requires  the  use  of  lifting-surface 
theory  to  calculate  the  aerodynamic  center  on  the  strut. 

Fourth,  structural  failure  can  occur  in  either  bending  or  torsion,  depending  on  the 
particular  stresses  involved.  This  happens  because  the  flow  produces  both  a  lift  and  a 
moment,  each  of  which  is  influenced  by  variation  in  angle  of  attack.  Prediction  of  actual 
failure  conditions  requires  a  detailed  stress  calculation  for  the  structure  as  a  function  of 
speed. 


FLAP  REVERSAL 

Flap  reversal,  a  form  of  control  surface  reversal,4  is  the  reversal  in  direction  of  the  in¬ 
crement  in  lift  induced  by  a  flap  deflection.  The  reversal  occurs  when  the  effect  on  lift  of 
foil  elastic  deflection  cancels  the  effect  of  a  flap  deflection.  Specifically,  when  a  trailing- 
edge  flap  is  deflected  downward,  the  resulting  hydrodynamic  twisting  moment  tends  to 
reduce  the  foil  angle  of  attack.  Since  the  torsional  rigidity  of  the  foil  is  independent  of 
speed  but  the  twisting  moment  increases  with  the  square  of  the  speed,  the  net  lift  increment 
will  be  zero  at  some  critical  speed,  namely  the  flap  reversal  speed.  This  condition  obviously 
would  limit  controllability  of  the  craft.  Both  the  absence  of  flap  reversal  on  full-scale  craft 
and  preliminary  design  calculations  for  the  PHM  craft  show  that  hydrofoils  have  flap  reversal 
speeds  well  above  the  subcavitating  speed  range.  Calculated  flap  reversal  speeds  are  con¬ 
sidered  sufficiently  accurate  for  subcavitating  hydrofoil  design.  Flap  reversal  on  cavitating 
foils  involves  substantially  different  loading,  different  foil  structures,  and  higher  speeds  and 
has  not  been  adequately  studied. 


RECOMMENDED  HYDROELASTIC  DESIGN  PROCEDURE 

For  desi^ ;  purposes,  a  T-foil  must  be  shown  to  be  free  from  both  flutter  and 
divergence  independently  of  each  other.  Both  modes  of  instability  must  be  investigated  for 
all  flow  conditions  because  the  speeds  at  which  failure  occurs  in  each  may  change  their 


4 

Scanlan,  R.H.  and  R.  Rosenbaum,  "Aircraft  Vibration  and  Flutn  ”  Dover  Publications,  Inc.,  New  York  (1968). 
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revive  positions  at  different  depths,  angles  of  attack,  etc.  If  a  model  is  used,  failure  may 
jccur  at  different  fractions  of  divergence  speed  on  mode!  and  prototype  because  of 
structural  differences. 


INITIAL  DESIGN 
St  ocavitating  Systems 

hi  view  of  the  large  number  of  stable  T-foils  in  operation,  initial  subcavitating  designs 
should  be  based  on  the  characteristics  of  existing  systems.  Such  systems  define  “design 
families”  composed  of  geometrically  similar  parts  and  similar  materials  which  will  also  be 
stable.  Only  the  size  must  be  specified,  according  to  the  loading  expected  from  craft  dis¬ 
placement,  maneuvering  requirements,  and  sea  conditions.  Existing  subcavitating  T-foils  are 
listed  in  order  of  increasing  chord  length  in  Table  1. 

If  existing  designs  are  not  suitable,  n  is  recommended  that  existing  configurations  be 
modified  in  the  light  of  known  parameter  trends. 


Cavitating  Systems 

All  hydrofoil  structures  wili  cavitate  at  spee  Is  above  50  knots.  No  design  information 
is  available  for  ensuring  stability  of  full-sized  T-foils  above  62  knots,  the  maximum  speed 
attained  by  the  DENISON.5  The  FRESH- 1  demonstration  T-foils  are  not  considered 
relevant  because  of  their  relatively  small  size,  as  indicated  in  Table  2.  The  design  of  cavi¬ 
tating  systems  requires  producing  a  structure  of  maximum  stability  under  subcavitating  con¬ 
ditions  and  carefully  avoiding  the  cavitation  patterns  which  are  known  to  be  destabilizing 
for  some  systems.  Initial  designs  for  operation  above  62  knots  must  be  regarded  as  highly 
tentative  until  their  stability  has  been  confirmed  by  model  experiments  or  operational 
experience. 


5Schaffer,  A.P.,  “Flutter  Flight  Testing  of  a  High  Speed  Hydrofoil  Craft-The  H  .  DENISON,”  Grumman  Aircraft 
Engineering  Corporation  Report  FT-HYD-la.l  (Apr  1963). 
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TABLE  I  -  EXISTING  SUBCAVITATING 
T-FOILS 


Craft 

Strut- 

Foil 

System 

Strut 

Chord 

in. 

PGH-2  (TUCI  MCAR  I 

Forward 

36 

PGH-1  (FLAGbTAt-F) 

Main 

40 

PGH-2  (TUCUMCARI) 

Aft 

45-55.6 

PCH  (HIGH  POINT) 

Forward 

51 

PGH-1  (FLAGSTAFF) 

Tail 

60 

AGFH  (PLAINVIEW) 

Tail 

72 

AG  EH  (PLAINVIEW) 

Main 

141 

TABLE  2  -  EXISTING  C'AVITATING  T-FOILS 


Craft 

Strut- 

Foil 

System 

Strut 

Chord 

in. 

FRESH-1 

Demonstration 

25.7 

Centerline 

FRESH-1 

Demonstration 

25.7 

Outboard 

DENISON 

Tail 

58.4-145 
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EVALUATION  OF  FIXED  DESIGN 


Three  techniques  can  be  used  to  estimate  the  flutter  and  divergence  stability  of  a 
proposed  T-foil  design.  It  is  helpful  to  an  analyst  to  use  as  many  of  the  techniques  as  possible 
in  studying  a  configuration. 


Model  Experiments 

Accurate  flutter  speed  estimates  can  be  obtained  from  experiments  with  a  flutter  model 
of  a  prototype.  Typical  subcavitating  prototype  structures  can  be  simulated  by  a  simplified, 
solid-sectior.  model  that  permits  useful  flutter  speed  estimates  to  be  obtained  at  moderate 
cost.  Ventilated  systems  may  require  a  model  which  is  slightly  more  complex.  Models  can 
be  built  at  considerably  greater  cost  for  studying  prototype  structures  in  detail  or  for 
studying  unusual  structures. 

The  use  of  models  for  determining  divergence  characteristics  would  not  be  productive 
unless  analytical  methods  are  invalidated  by  an  unusual  design  or  by  the  presence  of  signifi¬ 
cant  amounts  of  cavitation  or  ventilation. 

Operational  experience  of  the  stability  of  one  T-foil  in  a  design  family  may  be  regarded 
as  a  model  test  for  other  size  systems. 


Theoretical  Analysis 

Theoretical  predictions  of  flutter  speeds  of  T-foils  cannot  be  relied  on  since  adequate 
data  are  not  available  for  evaluation.  However,  theoretical  analysis  can  provide  valuable  in¬ 
formation  about  the  probable  flutter  mode  and  can  give  qualitatively  useful  guidance  on 
system  stability. 

Theoretical  divergence  predictions  for  subcavitating  flow  are  quite  accurate  and  suitable 
for  full  evaluation  of  system  stability.  Divergence  predictions  for  cavitating  flow  have  not 
been  evaluated. 


Full-Scale  Experiments 

Full-scale  flutter  experiments  provide  full  assurance  of  system  stability.  The  intent  is 
to  avoid  structural  failure  by  detecting  decreasing  damping  before  instability  occurs.  This 
may  be  fersible  for  torsional  flutter  but  may  not  be  for  bending  flutter. 
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Full-scale  flutter  testing  techniques  are  not  adequately  developed  at  the  present  time. 
Stability  measurements  require  a  method  of  structural  excitation  in  order  to  measure  damping. 
None  of  three  demonstrated  methods  appears  completely  satisfactory.  Impulsive  excitation 
by  log  impact,  achieved  randomly  and  unintentionally  on  the  AGEH  main  T-foil,  produced 
excellent  torsional  decay  signals  but  in  an  uncontrolled  manner.  The  DENISON  was  reported 
to  have  achieved  sufficient  excitation  by  maneuvering  and  by  crossing  a  wake.5  These 
techniques  may  be  inaccurate  because  the  forms  of  the  excitation  are  needed  for  proper 
analysis  and  these  are  not  known.  For  maximum  accuracy,  an  excitation  system  should  be 
attached  to  the  strut,  pod,  or  foil  or  else  a  known  impulsive  load  should  be  applied  ex¬ 
ternally.  Evaluation  of  flutter  stability  by  full-scale  testing  would  therefore  entail  development 
of  an  excitation  technique. 


SCALING  LAWS 

Scaling  laws  or  laws  of  similarity  must  be  applied  to  extrapolate  information  from 
existing  designs  or  model  studies  to  a  new  design.  It  is  possible  to  formulate  similarity 
conditions  for  two  systems  that  are  expected  to  be  kinematically  similar.  However,  it  is  not 
now  possible  to  state  any  general  stability  boundaries  in  quantitative  form  without  a  great 
deal  more  experimentation.  For  this  reason,  it  is  hoped  that  the  following  relationships  will 
be  useful  in  the  interim. 


SYSTEM  PARAMETERS 

System  parameters  are  requirtd  both  to  describe  a  given  T-foil  system  and  iu  provide  a 
basis  for  deriving  scaling  laws.  The  hydroelastic  analysis  of  T-foils  is  carried  out  by  using 
virtually  the  same  parameters  as  for  aeroelastic  anaiysis  of  airplane  wings6  but  with 
additional  parameters  for  free-surface  effects  and  cavitation.  Specifically,  hydrofoil  struts 
and  foils  can  be  represented  as  beams  by  using  straight  elastic  axes  determined  by  local  shear 
center  locations.  Structural  and  hydrodynamic  properties  can  be  given  as  spanwise  distri¬ 
butions  of  chordwise  sectional  properties.  Chordwise  sections  of  swept  structures  can  be 
taken  either  streamwise  or  in  a  direction  normal  to  the  elastic  axis;  axis-normal  sections  are 
used  in  this  report. 


6Bisplinghoff,  R.T.  et  al.,  “Aeroelasticity,”  Addison-Wesley  Publishing  Company,  Inc.,  Cambridge,  Massachusetts  (1955). 
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The  parameters  used  to  describe  a  T-foil  are  listed  in  the  appendix. 


NONDIMENSIONAL  SCALING  RATIOS 

T-foils  which  are  different  in  size,  material,  or  other  characteristics  will  have  equal 
nondimensional  flutter  and  divergence  boundaries  when  all  nondimensional  parameter  ratios 
are  the  same.  Under  these  conditions,  the  flutter  speeds  and  frequencies  obtained  for  one  of 
the  systems  can  be  used  to  determine  flutter  speeds  and  frequencies  for  the  other  systems. 
Divergence  speeds  can  be  similarly  scaled,  but  the  scaling  calculation  does  not  generally 
indicate  the  speed  of  structural  failure.  Nondimensional  flutter  and  divergence  parameter 
ratios  for  subcavitating  and  cavitating  T-foils  are  given  in  the  appendix.  The  scaling 
relationships  involving  system  size  and  stiffness  are  included  in  the  discussions  of  parameter 
trends. 


FLUTTER  OF  SUBCAVITATING  STRUT  SYSTEMS 
KNOWN  PARAMETER  TRENDS 

The  following  discussion  applies  primarily  to  subcavitating  systems  with  streamlined, 
airfoil-like  profiles.  The  parameters  of  some  cavitating  systems,  notably  those  with  blunt- 
based  profiles  which  are  naturally  ventilated  aft  of  the  section,  have  similar  trends. 
Differences  in  flutter  characteristics  of  subcavitating  and  cavitating  systems  are  discussed  in 
another  section  of  this  report. 


Flutter  Mode 

Strut  systems  undergo  both  bending  flutter  and  torsional  flutter.  Bending  flutter  occurs 
when  the  second  bending  frequency  is  lower  than  the  first  torsion  frequency  in  water. 
Torsional  flutter  occurs  when  the  reverse  is  true. 


Inertial  Parameters 


The  inertia  characteristics  of  a  T-foil  can  be  expressed  in  terms  of  generalized  mass  ratio 
li.  The  dependence  of  flutter  speed  on  n  is  shown  in  Figure  5.  The  most  stable  inertial 
configuration  for  a  strut  system  occurs  on  the  boundary  between  bending  flutter  and 
torsional  flutter. 


Bending  Flutter  Region.  As  shown  by  the  data  in  Figure  6,  flutter  speed  in  the  bending 
flutter  region  increases  monotonically  as  increases.1 ,7_1°  Therefore,  flutter  speedl 

increases  as  the  mass  of  the  strut  m  increases  and  as  the  submergence  C  decreases.  It  is 
further  known1  that  flutter  speed  increases  as  the  sweep  parameter  k  increases,  so  that  in¬ 
creasing  sweep  angle  A  or  decreasing  strut  length  L  will  stabilize  a  strut.  However,  bending 
flutter  probably  does  not  occur  when  struts  are  unswept,  implying  that  the  dependence  of 
flutter  speed  on  A  is  reversed  at  low  sweep  angles.  Experimental  evidence  for  these 
relationships  has  been  obtained  almost  entirely  from  struts  with  blunt-based  profiles. 


Torsional  Flutter  Region.  For  most  practical  hydrofoils  in  the  torsional  flutter  region, 
flutter  speed  decreases  rapidly  as  Mtorrion  increases.1  A  minimum  flutter  speed  wil'  always 
occur  if  submergence  is  decreased.  This  relationship  is  demonstrated  in  Figure  5  and  shown 
experimentally1 ,8,10_  12  in  Figure  7.  Therefore,  the  effects  on  flutter  speed  of  moment  of 
inertia,  pod  size,  and  other  parameters  in  Mtorsion  are  opposite  on  either  side  of  the  minimum 
flutter  speed.  It  is  expected,  however,  that  nearly  all  practical  configurations  will  be  to  the 
left  of  the  minimum  in  the  flutter  speed  curve. 


' 1 

Hilbome,  D.V.,  “The  Hydroelastic  Stability  of  Struts,”  Admiralty  Research  Laboratory  (Great  Britain)  Report  ARL / 
Rl/G/HY/5/3  (1958). 

8Huang,  T.T.,  “Experimental  Study  of  a  Low  Modulus  Flutter  Model  for  Strut-Foil-Pod  Configurations,” 

Hydronautics,  Inc.  Technical  Report  459-2  (Jul  1967). 

0  _ 

Squires,  C.E.,  Jr.,  “Hydrofoil  Flutter,  Small  Sweep  Angle  Investigation-Final  Report,”  Grumman  Aircraft  Engineering 
Corporation  Report  DA  Nonr-3989.3  (Nov  1963). 

'^Baird,  E.F.  et  al.,  “Investigation  of  Hydrofoil  Flutter-Final  Report,”  Grumman  Aircraft  Engineering  Corporation 
Report  DA  10  480-3  (Feb  1962). 

1  'Abramson,  H.N.  and  G.E.  Ransleben,  Jr.,  “An  Experimental  Investigation  of  Flutter  of  a  Fully  Submerged  Subcavitating 
Hydrofoil,"  J.  Aircraft,  Vol.  2,  No.  5,  pp.  439-442  (Sep-Oct  1965). 

i  2 

Besch,  P.K.  and  Liu,  Y.-N.,  “Flutter  and  Divergence  Characteristics  of  Four  Low  Mass  Ratio  Hydrofoils,”  NSRDC 
Report  3410  (Jan  1971). 
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(KNOTS) 


Figure  6  -  Experimental  Flutter  Speed  Uf  as  a  Function  of  Bending 
Mass  Ratio  Mbending  for  Solid-Section,  Blunt-Based  Struts 
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3  4 

^TORSION 


Figure  7  -  Experimental  Flutter  Speed  Uf  as  a  Function  of  Torsional 


Mass  Ratio  nt 


for  Torsion-Type  Struts 


Model  data111,12  for  both  simple  struts  and  for  a  strut  with  an  attached  pod  exhibited 
an  increase  in  flutter  speed  as  model  inertia  was  decreased;  see  Tables  3-6.  When  the  density 
of  the  smaller  strut  family  was  reduced,  the  flutter  speed  increased  from  24.7  knots  to  a 
value  above  36  knots  (Tables  3  and  4).  These  struts  remained  in  the  torsional  flutter  region 
at  low  densities  because  their  elastic  properties  were  concentrated  in  a  spar  of  unusually  low 
torsional  flexibility.  A  strut  with  a  large,  heavy  pod,  described  in  Tables  S  and  6,  exhibited 
an  increase  in  flutter  speed  when  pod  mass  was  decreased. 

It  is  useful  to  note  that  a  minimum  also  occurred  in  the  value  of  reduced  flutter  speed 
2Uf/ca>a  as  a  function  of  Minion  •  The  data  in  Figure  7  have  been  replotted  in  Figure  8  to 
illustrate  this  dependence.  Airfoil  data,13  also  plotted  in  Figure  8,  showed  similar  behavior. 
The  minimum  value  of  reduced  flutter  speed  for  hydrofoils  occurred  between  mass  ratios  of 
2.0  and  3.0.  The  lowest  value  of  reduced  flutter  speed  yet  obtained  for  streamlined  or 
blunt-based  struts  has  been  0.757.  A  rough  estimate  of  the  minimum  flutter  speed  for 
struts  with  elastic  properties  similar  tc  those  already  tested  can  be  made  from  the  relation 
2Uf/ccoa  >0.8  where  u>a  is  the  frequency  of  the  first  torsional  vibration  mode  in  air. 

X  Actual  flutter  speeds  can  be  many  times  larger  than  this  minimum. 


Elastic  Parameters 

The  elastic  characteristics  of  a  T-foil,  which  include  bending  stiffness  El,  torsional 
stiffness  GJ,  and  elastic  axis  location  xet,  affect  the  value  of  n  at  which  bending  flutter  and 
torsional  flutter  boundaries  intersect  as  well  as  the  magnitude  of  flutter  speed  in  each  region. 
Changes  which  move  the  second  and  third  vibration  mode  frequencies  closer  together  shift 
the  system  toward  the  boundary  between  bending  and  torsional  flutter. 


Size  and  Stiffness 

The  flutter  speed  of  a  given  configuration,  referred  to  as  the  model,  can  be  used  to 
determine  the  flutter  speed  of  a  different  configuration,  referred  to  as  the  prototype,  if  all 
dimensionless  parameters  given  in  the  appendix  are  equal  for  the  two  configurations.  Both 


13Woolston,  D.S.  and  G.E.  Castile,  "Some  Effects  of  Variations  in  Several  Parameters  Including  Fluid  Density  on  the 
Rutter  Speed  of  Light  Uniform  Cantilever  Wings,”  NACA  TN  2558  (1951). 
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TABLE  3  -  DIMENSIONED  PARAMETERS  FOR  SUBCAVITATING  STRUT 
FLUTTER  MODELS  IN  FRESH  WATER 


Strut 

Abramson  and 
Ransleben  Data11 

Besch  and  Liu  Data12 

Model  1 

Model  2 

Model  4 

L,  in. 

30.0 

15.0 

15.0 

15.0 

c,  in. 

12.0 

6.0 

6.0 

6.0 

A,  deg 

0.0 

0.0 

0.0 

0.0 

Xea'  in- 

3.0 

1.5 

1.5 

1.5 

V  ia 

6.14 

3.07 

3.07 

3.07 

m,  Ib/in. 

4.04 

0.983 

0.465 

0.206 

1  lb-in. 

my 

74.5 

4.49 

2.10 

0.938 

El,  lb-in.2 

3.40  x  106 

5.69  x  104 

5.53  x  104 

5.67  x  104 

GJ,  lb-in.2 

9.73  x  105 

1.92  x  104 

1.64  x  104 

1.56  x  104 

Profile: 

NACA  16-012 

NACA  16-012 

NACA  16-012 

NACA  16-012 

Flutter  Condition 

C.  in. 

30.0* 

15.0* 

15.0* 

15.0* 

Uf,  knots 

48.1 

24.7 

>36** 

>30 

ff,  Hz 

17.5 

20.5 

— 

— 

Mode  Shape 

Torsion 

Torsion 

— 

— 

Reflecting  plate  at  root. 

Structural  failure  due  to  approaching  divergence  occurred  at  36  knots. 


TABLE  4  -  NONDIMENSIONAL  PARAMETERS  FOR 
SU3CAVITATING  STRUT  FLUTTER  MODELS 


Strut 

Abramson  and 

Besch  and  Liu  Data 

,12 

Ranslebsn  Data11 

Model  1 

Model  2 

Model  4 

L/c 

2.5 

2.5 

2.5 

2.5 

AR 

5.0 

5.0 

5.0 

5.0 

K 

0.0 

0.0 

0.0 

0.0 

X.  /c 

0.25 

0.25 

0.25 

0.25 

fta 

0.51 

0.51 

0.51 

0.51 

4m/flp  c2 

0.99 

0.963 

0.455 

0.202 

iB^y/irpc4 

0.507 

0.489 

0.229 

0.102 

Flutter  Condition 

AR 

5.0 

5.0 

5.0 

5.0 

knf 

0.676 

0.772 

— 

— 

4  El/irp  c2  wf2  L4 

3.27  x  10-2 

2.31  x  10-1 

— 

— 

16  GJ/irp  L2  CJf2  c4 

2.34  x  10'1 

2.16  x  10_1 

— 

— 

Rnf 

6.7  x  106 

1.7  x  106 

— 

— 

Fnf 

oo 

OO 

oo 

oo 

^torsion 

1.35 

1.30 

0.610 

0.272 

^bending 

0.99 

0.963 

0.455 

0.202 

?s  (Mode  2) 

0.02 

0.02 

0.03 

0.05 

TABLE  5  -  DIMENSIONED  FLUTTER  PARAMETERS 
FOR  SUBCAVITATING  STRUT-POD 
FLUTTER  MODELS 

(From  Besch  and  Liu*) 


Strut 

1 - 

L,  in. 

58.4 

c,  in. 

12.0 

A,  deg 

15.0 

xea'  in‘ 

4.5 

Xcg'  *n‘ 

6.2 

m,  Ib/in. 

0.588 

W  lb-in’ 

8.21 

El,  lb-in.2 

7.58  x  105 

GJ,  lb-in.2 

10.4  x  105 

Profile 

NACA  16-005 

Pod 

Config.  A. 

Config.  B 

Length,  in. 

44.0 

44.0 

Max.  Diameter,  in. 

6.0 

6.0 

xn  .  in. 
nose 

16.2 

16.2 

X  •  in. 

22.0 

eg 

M,  lb 

140.0 

77.3 

ly,  lb-in.2 

1.57  x  104 

7.07  x  103 

Flutter  Condition 

C,  in. 

30.4 

30.4 

Uf,  knots 

6.6 

8.4 

ff,  Hz 

2.9 

3.5 

Mode  Shape 

Torsion 

T  orsion 

TABLE  6  -  NONDIMENSIONAL  FLUTTER 
PARAMETERS  FOR  SUBCAVITATING 
STRUT-POD  FLUTTER  MODELS 

(From  Besch  and  Liu*) 


Strut 

L/c 

4.87 

AR  (to  root) 

4.54 

K 

0.0551 

*ea/c 

0.375 

Xeg/C 

0.52 

4m/ffp  c2 

0.144 

16 

0.0559 

Pod 

Config.  A 

Config.  B 

Fineness  Ratio 

7.33:1 

7.33:1 

Diameter/c 

0.50 

0.50 

Xnt*/c 

1.35 

1.35 

x/c 

eg 

1.83 

1.83 

M  7M 

pod  strut 

4.08 

2.25 

*y,pod^*y, strut 

32.7 

14.7 

Flutter  Condition 

AR  (strut  wetted  area) 

2.36 

2.36 

knf 

0.846 

0.802 

4  El/trp  c2  uf2  L4 

0.0186 

0.0127 

16  GJ/np  L2  wf2  c4 

2.41 

1.66 

Rnf 

9.5  x  105 

1.2  x  106 

Fnf 

1.9 

2.5 

fs  (Mode  2) 

0.02 

0.02 

^torsion 

4.75 

2.21 

^bending 

1.15 

0.734 

2U,/ccj. 


AIRFOIL  MODELS 

(WOOLSTON  AND  CASTILE13) 

HYDROFOIL  MODELS 

EXPERIMENT 

MODEL 

EXPERIMENT  MODEL 

Q 

39-43-1 

• 

BESCH  AND  LIU12 

~  Ci 

39-42-3 

A 

HUANG8 

Q 

27-3B-4 

♦ 

BAIRD10 

0 

27-38-3 

• 

ABRAMSON  AND  RANSLEBEN 

-  V 

27-38-2 

• 

BESCH  AND  LIU17 

□ 

27-39-1 

i  ° 

17-32-4 

L  * 

17-32-3 

r  o 

17-32-2 

r~ 

[•  » 


X  |BB 

—  BLE 
(CAVITATION) 


(LIMITED  BY 
DIVERGENCE) 


LINES  CONNECT  DATA  OBTAINED  FROM 
1  THE  SAME  ELASTIC  CONFIGURATION 

NACA  *=  NACA  16-005 

BB  »  BLUNT  BASE  WITH  FULLY  WETTED  CHORD 
BLE  *  BLUNT  LEADING  EDGE 


^TORSION 

Figure  8  -  Reduced  Flutter  Speed  2Uf/cooa  as  a  Function  of  Torsional 
Mass  Ratio  Mtorsion  for  Torsion-Type  Struts 
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sire  and  bending  or  torsional  stiffness  can  be  chosen  arbitrarily  for  a  subcavitating  prototype. 
Blunt-based  end  other  ventilated  configurations  are  discussed  further  in  the  section  on  flutter 
of  cavitating  strut  systems.  The  flutter  speeds  of  the  model  and  prototype  will  be  related  as 
follows: 


prototype 

c2 

c  model 

/^prototype  ^model 

C2 

''model 

/  ^  prototype  ^model 

model 

C2  \ 

v  prototype  f 

f  ^  model  ^prototype 

or  « 

c2  V 

prototype! 

1  ^model  ^prototype 

where  values  of  c  and  El  or  GJ  are  taken  from  equivalent  locations  in  the  two  systems. 

Effects  of  Foils 

T-foils  have  been  studied  much  less  than  struts  and  strut-pod  systems.  It  is  believed  that 
the  effect  of  foils  can  be  explained  within  the  framework  of  bending  and  torsional  flutter. 

An  attached  foil  radically  changes  the  inertial  characteristics  of  the  strut-fluid  system.  System 
elasticity  is  less  suongly  affected.  Therefore  a  foil  could  change  the  mass  ratio  of  the  system 
and  yet  have  little  effect  on  the  flutter  boundary  (shown  in  Figure  5)  that  is  characteristic 
of  the  strut. 

At  present,  T-foil  system  characteristics  can  be  estimated  only  qualitatively  fiom  system 
mode  shapes  and  frequencies.  Strong  coupling  of  mode  shapes  places  T-foils  generally  in  a 
transition  region  between  bending  and  torsional  flutter.  Moderate-to-high  aspect  ratio  foils 
move  the  second  bending  frequency  of  a  torsion-type  strut  close  to  the  first  torsion  frequency. 
This  result  suggests  that  the  foil  will  raise  the  flutter  speed  of  a  torsion-type  strut.  The 
flutter  speed  of  a  bending-type  strut  might  be  lowered  by  a  foil.  Both  increases  and  decreases 
in  flutter  speed  were  observed  by  Huang8  when  a  large  foil  was  attached  to  a  strut  with  a 
blunt-based  profile.  These  results  are  discussed  further  in  the  section  dealing  with  cavitating 
systems. 


IMPLICATIONS  FOR  PROTOTYPE  SYSTEMS 

The  principal  implication  of  available  results  is  that  within  the  size  range  of  existing 
craft,  several  T-foils  define  structural  design  families  which  are  flutter  free  throughout  the 
subcavitating  speed  range  or  up  to  SO  knots  at  small  angles  of  attack.  This  conclusion  is 
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based  on  observed  stability  characteristics  of  both  model  and  full-scale  T-foils.  Scaling  laws 
required  for  extending  the  observed  characteristics  are  supported  by  a  model  experiment 
which  successfully  reproduced  scaling  predictions  for  a  system  of  different  size.  These 
results  are  discussed  below. 


Validity  of  Kinematic  Scaling 

The  flutter  testing  of  two  nondimensionally  similar  simple  strut  models  has  prov  ded 
evidence  for  the  validity  of  kinematic  scaling  procedures  for  determining  flutter  bounda'ies. 
The  two  models111214  were  rectangular,  torsion-type  struts  o' aspect  ratio  5.0.  Mojel 
parameters  were  those  for  the  30-inch  model  and  Model  1  given  in  Tables  3  and  4.  The 

scaled  flutter  speed  of  the  larger  strut  agrees  with  its  experimental  flutter  speed  within  1  per¬ 
cent  when  El  is  used  for  scaling,  and  within  9  percent  when  GJ  is  used.  Exact  agreement  is 

not  expected  because  GJ/EI  and  4m/7rp  c2  were  not  exactly  equal  for  the  two  stmts. 


AGEH  Main  Design  Family 

The  AGEH  main  T-foil,  described  in  Tables  7  and  8,  is  now  established  to  be  free  from 
flutter  throughout  the  subcavitating  speed  range.  It  therefore  serves  as  a  parent  for  a  family 
of  T-foils  of  different  size  which  are  also  flutter  free.  The  characteristics  of  this  design  family 
will  be  described  following  the  discussion  of  experimental  results. 

Flutter  stability  of  the  AGEH  main  T-foil  has  been  established  by  both  full-scale  data 
and  1/8-scale  model  data.  The  full-scale  data  consisted  of  stable  craft  operation  through 
47  knots  plus  damping  values  at  two  intermediate  speeds  measured  from  oscillations  excited 
by  impacts  with  logs.*  The  measured  damping  values,  shown  in  Figure  9  at  35  and 
45.5  knots,  were  sufficiently  high  that  flutter  should  not  occur  in  the  torsional  (second) 
mode  until  above  50  knots.  Torsional  flutter  is  expected  to  precede  bending  flutter  in  this 
torsion-type  strut  system. 

A  1/8-scale  flutter  model15  of  the  main  foil  system,  constructed  as  part  of  the  con¬ 
tract  requirements,  also  provided  evidence  of  flutter  stability.  Parameters  of  the  model  are 


14 

Ranslebcn,  G.E.,  Jr.,  "Experimental  Delerminalion  of  Variation  of  Hydrofoil  Flutter  Speed  wilh  Mass  Ratio, '  South¬ 
west  Research  Institute,  Contract  N00014-69-C-0219  (Apr  1970). 

15Ranslcbcn,  G.E.,  Jr„  "Description  of  a  One-Eighth  Scale  Flutter  Mocel  of  the  Main  Hydrofoil  Assembly  of  the  AG(EH) 
800  Hydrofoil  Research  Ship,”  Southwest  Research  Institute  Project  02-1  *60  (Dec  1966). 

*AC.HH  Full-Scale  Trial  71HA2 1 1-22  (24  Aug  1971,  Voyage  66,  Tape  AT6bB,  Time  16:22:27)  aid  Trial  72HA-A241G-31 
(31  May  1972,  Tape  180B,  Voyage  117,  Time  10:39:30). 
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TABLE  7  -  DIMENSIONED  FLUTTER  PARAMETERS  FOR  FULL-SCALE 
SUBCAVITATING  T-FOILS 


Strut 

PCH 

Mod  1 

PHM 

AGEH 

Tail 

AGEH 

Main 

L,  in. 

129.6 

220.2 

278.1 

328.6 

!  c,  in. 

51.0 

66.5 

72.0 

141 

A.  deg 

0.0 

6.3 

0.0 

0.0 

xw  sP*n*’  in- 

14.8 

27-33 

43.2 

45 

x  (3/4  span),  in. 
eg 

23.9 

28.3 

36 

66 

m  (3/4  span),  Ib/in. 

25 

24 

19 

60 

1  (3/4  span),  lb-in, 

5.2  x  103 

1.3  x  104 

9.7  x  It^ 

1.0  x  105 

El  (3/4  span),  lb-in.2 

7x10® 

2.0  x  1010 

1.5  v  1u10 

1.1  x  1011 

GJ  (3/4  span),  lb-in,2 

9.5  x  109 

2.9  x  10’° 

1.5  x  1010 

1.1  x  1011 

Profile  (3/4  span) 

NACA  16-012 

NACA  16-012 

NACA  16-011 

NACA  16-012 

Pod 

Length,  in. 

oV 

157.5 

156 

321.7 

Max.  Diameter,  in. 

13.5 

25 

30 

49 

xn  .  in. 

nose 

24.0 

40 

41.2 

107.1 

c 

8 

X 

186 

35.4 

40 

97 

M,  lb 

226 

2,351* 

875 

13,705 

ly,  lb-in.2 

4  6  x  104 

7.4  x  106* 

1.98  x  106 

1.25  x  108 

Foil 

Area,  ft2 

65.6 

140.9 

SO.O 

225.0 

c  in. 

root 

63.0 

93.3 

82.4 

159.8 

Tip  Chord,  in. 

15,7 

28.0 

24.7 

47.9 

Full  Span,  in. 

240 

335 

161 

312 

A.  deg 

15.0 

11.04 

35.2 

35.2 

xfoil-  in- 

12.0 

13.2 

6.79 

51.7 

M.  lb 

4200 

7201 

2182 

8884 

xcg'  in- 

25  3 

37.8 

48 

58 

ly.  Ib-in.2 

1.16  x  107 

2  62  x  106 

4.94  x  107** 

Profile 

NACA  16-309 

NACA  16- 

NACA  16X08 

NACA  16X08 

a  =  1.0 

206.5 

x  =  0.425 

x  =  0.39 

a  -  1.0 

a  =  1.0 

a  ’  1.0 

Flutter  Condition 

8.  in. 

Not  available 

Not  available 

— 

96 

U(,  knots 

Above  50 

Above  50 

f(.  Hz 

— 

— 

Mode  Shape 

.lot  available 

Not  available 

— 

— 

*  Pod  inertie  includes  section  of  strut  enclosed  within  nod  and  an  estimated  558  lb  of  enclosed  water 
(pod  is  flooded!. 


•  • 

About  vertical  axis;  strut  elastic  axis  has  10-deg  dihedral. 
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TABLE  8  -  NONDIMENSIONAL  FLUTTER  PARAMETERS  FOR  FULL-SCALE 

SUBCAVITATING  T-FOILS 


AR  (to  root) 


xea  ^/4  span)/c 
(3/4  span)/c 
4m  (3/4  span)/trp  c2 
16  lmy  (3/4  span)/7Tp  c4 


Fineness  Ratio 
Diameter/c 

Xno*/c 
%,pod/c 
Mpod/Mstru, 
^y.pod^y. strut 


C  /c 
root' 


A,  deg 


Mfoi,/Mstrut 
*y,foH^y,ftrut 
Flutter  Condition 


AR  (strut  wetted  area) 


4  El  (3/4  span)/7rp  c2  wf2  L4 
16  GJ  (3/4  span)/irp  L2  OJf2  c4 


AGEH 

Tail 


0.40-0.50 


0.0941 


0.0497 


AGEH 

Main 


0.0348 


No*  available  Not  available  Not  available  Not  avai'able 


Not  available 

move  8.1  x  107 

i  above  4.3 

' 

Not  available  Not  available  Not  available  0.09  (Mode  2) 


$ 


SPEED  U  (KNOTS) 


Figure  9  -  Hydroelastic  Mode  Characteristics  for  AGEH  Main  T-Foil 
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given  in  Tables  9  and  10.  This  model  (shown  in  Figure  1)  was  constructed  of  single  spars 
with  independent  skin  sections  and  involved  considerable  expense.  The  scale  model  was 
designed  by  selecting  the  desired  length  and  velocity  ratios  between  prototype  and  model,  thus 
determining  model  stiffnesses.  In  this  case,  the  velocity  ratio  was  chosen  to  produce  the  same 
value  of  Froude  number  Fn  for  both  systems.  This  condition  was  not  required,  however,  at 
the  relatively  high  values  of  Fn  associated  with  the  system. 

Model  vibration  characteristics  were  in  fairly  good  agreement  with  full-scale  character¬ 
istics  but  did  exhibit  some  discrepancies.16  The  model  was  relatively  more  flexible  in 
torsion  than  the  prototype  and  had  2  percent  less  structural  damping  in  the  second,  or 
unstable  mode.  When  the  model  frequency  was  appropriately  sealed,  the  greater  flexif '  ty 
of  the  model  was  revealed  by  the  fact  that  its  torsional  frequency  (at  29-percent  sub¬ 
mergence)  was  lower  than  that  of  the  prototype  (at  80-percent  submergence). 

The  structural  damping  in  the  second  mode  of  both  model  and  prototype,  7  and  9  per¬ 
cent,  respectively,  was  unusually  high  for  winglike  structures.  Such  high  damping  raises  the 
flutter  speed  of  the  second  mode  substantially.  It  is  encouraging  that  both  model  and  proto¬ 
type  developed  high  Mode  2  damping  despite  differences  in  construction  because  it  suggests 
that  structural  damping  is  to  some  extent  simulated  even  though  it  is  usually  not  a  design 
parameter. 

Figure  9  shows  model  frequency  data  obtained  from  mechanical  impedance  measure¬ 
ments  made  at  NSRDC  during  flutter  testing  (unpublished  data  from  D.S.  Cieslowski).  Fre¬ 
quencies  have  been  converted  to  full-scale  values.  Damping  values  obtained  from  impedance 
measurements  showed  excessive  scatter  and  are  not  given.  The  damping  ratio  at  zero  speed 
was  obtained  by  line  cut  excitation. 

At  33  knots,  the  model  encountered  a  flutter  instability  at  a  frequency  equal  to  that  of 
the  second  mode  at  zero  speed.  The  strut  was  29-percent  submerged,  corresponding  to 
design  flying  height.  The  instability  damaged  the  model  and  precluded  testing  at  other  sub¬ 
mergence  depths.  Previous  results1  indicate  that  minimum  torsional  flutter  speeds  of  strut- 
pod  models  occurred  at  50-percent  submergence  but  that  one  T-foil  model  had  a  minimum 
flutter  speed  at  100-percent  submergence.  It  is  concluded  that  the  minimum  flutter  speed  of 
the  AGEH  model  is  probably  lower  than  33  knots  and  that  it  would  occur  at  a  deeper  sub¬ 
mergence. 

This  model  result  correpsonds  to  a  scaled  flutter  speed  of  93  knots  for  the  full-scale 
AGEH  main  T-foil  system.  This  prediction  is  valid  only  for  subcavitating  flow  conditions 


6Peoples,  J.R.,  “Frequencies  and  Damping  of  Full  Scale  Hydrofoils  by  “Pluck  Test”  Methods,”  43rd  Shock  and 
Vibration  Bulletin,  Shock  and  Vibration  Information  Center,  Naval  Research  Laboratory,  Washington,  D.C.  (Jun  1973); 
also  indicated  by  D.S.  Cieslowski  of  NSRDC  in  a  private  communication. 
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TABLE  9  -  DIMENSIONED  FLUTTER  PARAMETERS  FOR 
1/8-SCALE  AGEH  MAIN  T-FOIL  FLUTTER  MODEL 


Strut 

Foil 

L,  in. 

41.1 

Area,  in.2 

506 

c,  in. 

17.6 

in* 

root 

20.0 

A,  deg 

0.0 

Tip  Chord,  in. 

6.0 

x  (3/4  span),  in. 

5.6 

Foil  Span,  in. 

39.0 

x^  (3/4  span),  in. 

8.3 

A,  deg 

35.2 

m  (3/4  span),  in. 

0.916 

XfoiC  in- 

6.46 

1  (3/4  span),  lb-in. 

my 

2.98 

M,  lb 

16.9 

El  (3/4  span),  lb-in.2 

3.28  x  106 

V  in- 

7.3 

GJ  (3/4  span),  lb-in.2 

3.28  x  106 

ly,  Ib-in.2 

1.47  x  103* 

Profile  (3/4  span) 

Pod 

NACA  13-012 

V  Ib-in.2 

Profile 

NACA  16X08 

x  =  0.390 
a  =  1.0 

Length,  in. 

40.2 

Flutter  Condition 

Max.  Diameter,  in. 

6.1 

2,  in. 

12.0 

Xnose'  in‘ 

13.4 

Uf,  knots 

33 

w in- 

12.1 

ff,  Hz 

'j  0.6 

M,  lb 
ly.  Ib-in.2 

V  lb-!"-2 

26.1 

3.73  x  103* 

Mode  Shape 

Torsion 

* 

About  vertical  axis;  strut  elastic  axis  has  10-deg  dihedral. 


TABLE  10  -  NONDIMENSIONAL  FLUTTER  PARAMETERS  FOR 
1 /8-SCALE  AGEH  MAIN  T-FOIL  FLUTTER  MODEL 


Strut 

Foil 

L/c 

2.34 

c  „./c 

root 

1.K 

AR  (to  root) 

2.3 

AR 

KTh 

K 

0.0 

A 

35.2 

xea  (3/4  span)/c 

0.32 

T 

0.30 

xcg  (3/4  span)/c 

0.47 

xfoil/c 

0.3S7 

4m  (3/4  span)/irp  C2 

0.104 

Wc 

0.41 

16  1  (3/4  span)/irp  c4 

my 

0.0348 

Mfo,/Mstrut 

0.42 

(y.foil^v, strut 

0.59 

Pod 

’p.foil'  (p.pod 

Fineness  Ratio 

6.6:1 

Diameter/c 

0.347 

Flutter  Condition 

x  /c 
nose 

0.761 

AR  (strut  wetted  area) 

0.68 

Vs 

0.688 

knf 

0.88 

Mpod/Mstru, 

0.65 

4  El  (3/4  span)/irp  c2  to,2  L4 

0.0114 

(y,poc/(y  .strut 

1.48 

16  GJ  (3/4  span)/7 rp  L2  cof2  c4 

0.247 

(^.pod^y, strut 

Rn. 

6.7  x  106 

Fnf 

8.1 

f. 

0.07 

and  a  strut  submergence  of  29  percent.  Subcavitating  flow  can  be  maintained  only  up  to 
approximately  50  knots  at  small  angles  of  attack;  the  effect  of  cavitation  on  flutter  of  stream¬ 
lined  profiles  has  not  been  studied.  Oil  the  basis  of  previous  results,1  it  appear;  unlikely  that 
flutter  speed  would  decrease  as  much  as  46  percent  as  submergence  increased  Furthermore, 
the  small  discrepancies  in  damping  and  torsional  stiffness  between  model  and  prototype  would 
raise  the  scaled  flutter  speed.  Therefore  the  AGEH  main  T-foil  is  considered  flutter  free  for 
all  subcavitating  flow  conditions. 

A  family  of  similarly  flutter-free  T-foil  designs  are  obtainable  from  this  result.  The 
family  consists  of  all  T-foils  constructed  of  components  geometrically  scaled  from  the  AGEH 
system  and  composed  of  the  same  (or  equal  density  and  moduli)  materials.  The  section 
moduli  of  such  T-foils  vary  as  the  fourth  power  of  any  linear  dimension.  Therefore,  if  the 
structural  damping  ratios  of  the  second  modes  are  equal,  the  scaled  systems  will  flutter  at 
the  same  speed  as  the  AGEH  system,  according  to  the  scaling  relationship  presented  earlier. 
This  should  be  confirmed  by  full-scale  vibration  tests  of  future  systems.  Froude  number  is 
not  preserved  under  these  conditions,  however,  and  it  is  expected  that  scaling  would  become 
inaccurate  for  extremely  large  systems  which  correspond  to  low  values  of  Froude  number. 
Torsional  stiffnesses  GJ  for  geometrically  scaled  systems  correspond  to  the  AGEH  main 
design  fir.uiy  shown  in  Figure  10.  A  more  efficient  structural  design  would  give  greater 
rigidity  for  the  same  weight  of  material  and  thus  be  more  stable. 


Other  Stable  Desinn  Families 

The  AGEH  main  T-foil  serves  as  a  parent  for  a  design  family  of  stable  T-foils  with 
relatively  heavy  pods  and  moderate  aspect  ratio  foils.  All  other  existing  T-foils  similarly 
define  different  stable  families  of  subcavitating  T-foils  by  virtue  of  having  successfully 
operated  at  high  subcavitating  speeds.  Parameters  for  three  such  T-foils  are  given  in  Tables  7 
and  8.  The  rule  of  geometric  and  material  similaiity  can  be  used  to  produce  equally  stable 
T-foils  up  tc  a  size  at  which  Froude  number  is  considered  too  low  to  be  disregarded.  The 
design  family  for  the  PCH  and  PHM  struts,  shown  in  Figure  10,  gives  further  indication  that 
maximum  stiffness  for  operational  T-foils  has  probably  been  reached.  Intermediate  stiffness 
families,  corresponding,  e.g.,  to  the  AGEH  tail  T-foil,  can  be  established  by  drawing  lines 
parallel  to  the  lines  given  in  Figure  10. 

The  stability  of  all  operational  systems  suggests  that  present  design  criteria  that  are  not 
based  on  hydroelastic  considerations  do  yield  T-foils  which  are  sufficiently  rigid  to  preclude 
flutter  and  divergence.  The  design  criterion  for  the  PCH  Mod-1  forward  strut  was  the 
ability  to  withstand  full  loading  on  one  side  of  a  half-broached  foil  or,  alternatively,  a  0.2-g 


Figure  10  -  Torsional  Stiffness  GJ  as  a  Function  of  Strut  Chord 
Length  for  Geometrically  Similar  T-Foil  Designs 
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strut  side  load.17  Thus  although  present  designs,  particularly  those  for  the  smaller  chord 
PCH  and  PHM,  might  be  overdesigned  for  hydioelastic  stability,  other  considerations  would 
prevent  weight-saving  optimization. 

It  should  be  observeo  that  as  size  increases,  the  prototype  stiffnesses,  shown  in 
Figure  10,  increase  more  slowly  than  geometrically  scaled  stiffnesses.  This  trend  is  due  to 
structural  differences  in  the  systems  represented.  These  systems  do  not  necessarily  represent 
the  only  design  alternatives  as  a  function  of  size.  However,  if  large  size  struts  are  con¬ 
strained  by  fabrication  requirements  or  other  considerations  to  continue  this  trend,  then  the 
flutter  speeds  of  large  members  of  a  design  family  will  decrease  as  size  increases.  Therefore, 
existing  design  families  cannot  be  assumed  to  be  hydroelastically  stable  and  suitable  for  use 
on  craft  larger  than  the  AGEH  until  a  study  is  made  of  the  engineering  aspects  of  large-scale 
strut-foil  structures,  including  the  effects  of  size  on  other  design  criteria. 


EXPERIMENTAL  MODELING  TECHNIQUES 

Flutter  models  are  extremely  valuable  for  confirming  the  stability  of  hydrofoil 
structures.  Because  of  the  present  inadequacy  of  flutter  theory,  for  example,  quantitative 
bending  flutter  boundaries  can  be  obtained  only  from  model  results.  If  there  is  any  question 
of  flutter  instability  after  a  comparison  with  existing  designs  and  the  performance  of  a 
modified  strip  analysis,  a  model  experiment  should  be  conducted  before  construction  is 
started. 

Difficulties  in  designing  successful  strut  models  have  been  experienced  in  the  past,  but 
these  difficulties  can  now  be  avoided.  It  is  proposed  to  use  a  solid  cross  section  model  which 
is  a  dose  but  not  exact  simulation  of  a  prototype  configuration.  Model  requirements  are 
discussed  below. 


Model  Cost 

Model  cost  depends  on  the  complexity  of  construction.  Semimonocoque  models  can 
simulate  all  prototype  elastic  and  inertia  characteristics  exactly,  but  at  considerable  expense. 
It  is  therefore  proposed  to  use  solid  cross  section  models,  which  are  relatively  inexpensive  to 
build. 


17“PCH-1  Mod  1  Detail  Design-Structural  Summary,”  The  Boeing  Company  (1969). 
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Solid  Cross  Section  Parameters 

Solid  cross  sections  of  streamlined  profile  can  be  used  to  adequately  simulate  the  elastic 
and  inertial  properties  of  most  prototype  struts  and  foils  under  subcavitating  flow  conditions. 
The  simulation  is  based  on  choosing  the  thickness-to-chord  ratio  of  the  model  so  that  the 
required  mass  and  moment  of  inertia  are  produced.  The  prototype  thickness  does  not  have 
to  be  simulated  exactly  so  long  as  cavitation  does  not  occur. 

The  remaining  parameters  that  must  be  simulated  are  inherently  well  modeled  by  a 
solid  cross  section,  namely,  center  of  gravity  location  xcg/c,  elastic  axis  location  xea/c.  and 
GJ/E1.  Characteristic  values  of  these  parameters  for  existing  prototype  craft  and  for  a  solid 
cross  section  are  given  in  Table  1 1.  Although  prototype  values  of  xea/c  va-y  widely  on 
either  side  of  the  solid  section  value,  the  extreme  values  shown  are  mislead  ng.  The  farthest 
aft  location  of  60-percent  chord,  found  on  the  AGEH  tail  strut,  is  expected  to  be  avoided  in 
future  designs  because  it  leads  to  unfavorable  divergence  characteristics.  The  farthest  forward 
location  of  20-percent  chord  is  a  local  characteristic  produced  by  terminating  the  root  of  the 
strut  in  a  kingpost  for  steering  purposes;  this  can  be  simulated  by  a  similar  suspension.  The 
elastic  axes  of  subcavitating,  semimonocoque  struts  and  foils  are  expected  to  fall  between 
40-  and  50-percent  chord.  It  appears  that  some  variation  in  GJ/EI  may  have  to  be  accepted 
for  some  structures. 

The  ratio  of  the  flutter  speed  of  the  prototype  to  that  of  the  model  can  be  determined 
by  using  the  scaling  principles  described  earlier  which  relate  stiffness  and  flutter  speed.  For  a 
given  prototype,  this  ratio  will  depend  on  the  model  material.  For  one  system  studied,  steel, 
titanium,  and  aluminum  produced  velocity  scale  ratios  of  4.0,  2.1,  and  1.5,  respectively.  Any 
of  these  materials  would  be  suitable  for  studying  flutter  characteristics  of  a  subcavitating 
hydrofoil  system  in  a  50-knot  test  facility. 


Model  Size 

Models  must  undergo  flutter  under  flow  conditions  which  produce  forces  similar  to 
those  experienced  by  prototype  systems.  The  similarity  oi  oscillatory  forces  for  airfoil 
systems  of  different  sizes  can  be  approximated  for  values  of  Reynolds  number  Rn  above 
4  x  105,  as  indicated  on  page  710  of  Bisplinghoff  et  al,6  Hydrofoil  prototypes  have  values 
of  Rn  far  above  this  value,  and  models  of  practical  size  readily  exceed  it.  For  example,  a 
6-inch  chord  model  at  14  knots  has  a  value  of  Rn  =  10* 
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TABLE  1 1  -  COMPARISON  OF  PARAMETERS 
FOR  SOLID  CROSS  SECTION  MODELS  AND 
FULL-SCALE  STRUTS  AND  FOILS 


Range  of 

Value  for 

Parameter 

Prototype 

NACA  16-Series 

Values 

Solid  Section 

Center  of  Gravity 
Location,  x^/c 

0.46-0.50 

0.48 

Elastic  Axis 
Location,  xga/c 

0.20-0.60 

0.48 

GJ/EI 

1.0-1.45 

1.26  (aluminum) 
1.49  (steel) 

4m  (3/4  span) 

0.1-0.34 

same* 

irp  c2 

• 

Using  the  following  model  materials  and  thickness-to- 
chord  ratios: 

Aluminum,  4.1  to  13.4  percent  thick, 

Titanium,  2.3  to  7.4  percent  thick, 

Steel,  1.4  to  4.6  percent  thick. 

39 


Experimental  Configurations 
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All  flutter  models  must  be  tested  with  sufficient  sweepback  to  preclude  divergence. 
Divergence  speeds  can  be  calculated  with  sufficient  accuracy  (as  discussed  elsewhere  in  this 
report)  to  make  this  possible.  A  sweep  angle  of  10  degrees  is  usually  sufficient.  Models  with 
unswept  struts  common  to  full-scale  craft  will  often  diverge  before  experiencing  flutter,  but 
the  prototype  might  achieve  speeds  closer  to  divergence  and  thereby  undergo  flutter. 
Extrapolation  to  smaller  sweep  angles  should  be  made  theoretically. 

Models  with  load-bearing  foils  must  operate  well  below  the  critical  buckling  load  of  the 
supporting  strut.  Vertical  loading  of  near-buckling  magnitude  has  decreased  the  speed  at 
which  predivergence  deflections  occurred.8  Furthermore,  load-bearing  foils  can  cause  bend¬ 
ing  deflections  of  a  strut  and  make  testing  impossible  even  though  divergence  or  buckling  is 
not  involved.  Conventionally  propoitioned  foils  can  readily  produce  excessive  vertical  load¬ 
ing  on  model  struts  because  the  model  struts  have  relatively  lower  rigidity  than  prototype 
struts  in  order  to  simulate  speeds  above  the  prototype  operating  range.  Therefore  it  is 
recommended  that  models  permit  control  of  vertical  lift. 

It  is  desirable  to  prepare  pod  and  foil  in  additional  configurations  as  well  as  in  that 
corresponding  to  the  prototype  so  that  flutter  can  be  obtained  at  lower  speeds  than  for  the 
configuration  of  interest.  These  will  provide  additional  quantitative  flutter  boundaries  to  be 
used  in  combination  with  the  theoretical  projection.  The  additional  results  are  particularly 
important  if  the  prototype  configuration  does  not  undergo  flutter  within  the  speed  range  of 
the  test  facility.  Failure  of  the  prototype  configuration  to  flutter  within  the  50-knot  speed 
range  of  the  high-speed  towing  basin  at  NSRDC18  would  normally  correspond  to  a  flutter 
speed  for  a  similarly  swept  full-scale  structure,  well  above  speeds  at  which  subcavitating  flow 
can  be  maintained.  Experimental  results,  shown  in  Figures  7  and  8  at  Utortion  =1.1  demon¬ 
strate  that  large  amounts  of  ventilation  are  destabilizing.  The  effect  of  moderate  amounts 
of  cavitation,  such  as  would  be  expected  during  operation  of  a  streamlined  strut  system  in 
the  vicinity  of  60  knots,  has  not  been  investigated. 

Experimental  configurations  should  include  several  sweep  angles  above  the  divergence 
boundary.  They  should  also  include  several  depths  because  full-scale  strut  systems  operate 
from  zero  to  1 00-percent  submergence  in  rough  water  and  cannot  be  considered  to  be 
limited  to  their  design  flying  height.  As  previously  discussed,  existing  data  show  that 
flutter  speed  changes  with  strut  submergence. 


1  fi 

Brownell,  W.F.  and  M.L.  Miller,  “Hydromechanics  Cavitation  Research  Far  ties  and  Techniques  in  Use  at  the  David 
Taylor  Model  Basin,”  Symposium  on  Cavitation  Research  Facilities  and  Techniques,  American  Society  of  Mechanical 
Engineers  (May  1964);  also  David  Taylor  Model  Basin  Report  1856  (Oct  1964). 
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During  experimental  work,  care  should  be  taken  to  accelerate  the  flutter  model  slowly 
and  to  increase  speed  in  small  increments.  When  damping  is  measured  frequently,  as  can  be 
done  by  using  the  rapid  sweep  excitation  method  described  below,  it  is  possible  to  investi¬ 
gate  several  speeds  on  a  given  test  run. 


Damping  Measurement  Techniques 

The  preferred  method  of  testing  flutter  models  consists  of  measuring,  insofar  as  possible, 
the  damping  of  the  three  lowest  hydroelastic  modes  throughout  the  speed  range  below 
flutter.  Knowledge  of  damping  permits  flutter-and  thus  the  possibility  of  structural  failure 
of  the  model -to  be  anticipated  and  assists  in  theoretical  and  practical  understanding  of  the 
instability. 

Three  damping  measurement  techniques  are  available,  but  all  are  subject  to  limitations. 
The  most  versatile  technique  is  rapid  sweep  excitation.  A  vibration  generator  is  used  to 
apply  a  sinusoidal  force  of  rapidly  varying  frequency  to  a  model.  The  frequency  is  swept 
from  slightly  below  to  slightly  above  a  known  resonance  frequency.  When  the  sweep  is 
complete,  the  force  is  cut  off.  Damping  values  are  obtained  either  from  the  decaying 
oscillation  following  excitation  or  from  a  phase  plane  plot  of  the  frequency  response 
function.19  When  the  former  method  of  analysis  is  used,  this  technique  is  not  suitable  for 
measuring  damping  above  approximately  1 0  percent  of  critical  damping. 

A  second,  similar  technique  is  known  as  line  cut  excitation.  A  line  is  attached  to  a 
model,  pulled  so  as  to  deflect  the  model,  and  then  cut.  This  method  can  produce  a  decay 
curve  for  a  mode  of  oscillation  having  up  to  1 0  percent  of  critical  damping. 

A  third  damping  measurement  technique  involves  measuring  the  mechanical  impedance 
of  a  model  at  a  resonance,  thus  giving  the  damping  directly.  A  vibration  generator  and 
response  pickup  are  required.  This  technique,  which  should  be  suitable  for  measuring  high 
values  of  damping,  requires  long  data  acquisition  times  and  lias  not  beer,  thoroughly  evaluated 
for  flutter  work. 

Damping  measurement  is  hampered  by  the  large  range  of  damping  values  encountered 
for  some  struts  as  well  as  by  the  rather  low  frequencies  of  some  hydroelastic  modes.  Hydro- 
elastic  modes  of  flutter  models  usually  having  damping  ratios  f  of  5  percent  or  less  at  zero 
speed.  This  damping  is  due  almost  entirely  to  internal  structural  damping.  The  mode  which 


19 

White,  R.G.,  “Use  of  Transient  Excitation  in  the  Dynamic  Analysis  of  Structures,”  The  Aeronautical  Journal  of  the 
Royal  Aeronautical  Society  Vol.  73,  pp.  1047-1050  (Dec  1969). 
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is  unstable  in  torsional  flutter  sometimes  increases  in  damping  only  slightly  before  decreasing 
to  zero  at  flutter  inception.  In  this  low  damping  region,  rapid  sweep  excitation  permits 
damping  measurements  to  be  made  readily.  The  unstable  mode  in  the  bending  flutter  and 
transitional  flutter  regions,  however,  is  very  highly  damped  at  intermediate  speeds  and  drops 
rapidly  to  zero  damping  in  a  short  speed  range  prior  to  flutter  inception.  This  mode  is 
difficult  to  detect  prior  to  flutter,  particularly  because  of  its  frequencies  which  are  usually 
below  1  Hz.  It  might  be  possible  to  determine  damping  values  for  the  bending  flutter  mode 
from  mechanical  impedance  measurements. 

The  three  excitation  methods  that  have  been  demonstrated  on  full-scale  systems-Iog 
impact,  maneuvering,  and  crossing  a  wake-are  similar  to  the  line  cut  technique  but  have  the 
additional  limitation  that  deflection  amplitude  is  uncontrolled  at  the  present  stage  of 
development. 


SUBCAVITATING  FLUTTER  THEORY 

The  accuracy  of  flutter  speed  predictions  for  struts  with  attached  foils  is,  in  general, 
nknown.  The  one  theory  evaluated  to  date  shows  only  fair  agreement  with  one  experi¬ 
mental  result.  Predictions  for  simple  struts  and  struts  with  pods  have  a  known  degree  of 
inaccuracy,  and  this  permits  calculations  to  be  corrected  to  give  good  flutter  speed  estimates 
for  torsion-type  struts. 


Flutter  Predictions  for  T-Foils 

It  is  impossible  to  evaluate  the  accuracy  of  flutter  predictions  for  subcavitating  T-foils 
because  insufficient  flutter  data  are  available  for  such  systems.  Only  one  flutter  point 
exists,  corresponding  to  the  AGEH  main  T-foil.  This  result  is  used  for  a  single-point 
evaluation  of  a  modified  strip  flutter  theory  in  the  following  discussion.  In  view  of  the  wide 
variations  in  flutter  speed  and  mode  shape  expected  for  practical  T-foils,  this  evaluation  is 
not  of  general  applicability. 

The  previously  mentioned  full-scale  and  model  data  for  the  AGEH  main  T-foils 
constitute  the  available  T-foil  flutter  data.  The  data  are  shown  in  terms  of  prototype  values 
in  Figure  9. 


Among  applicable  flutter  theories,  only  the  modified  strip  flutter  theory  has  been  widely 
studied1 ,20~ 22  and  can  be  considered  a  possible  candidate  for  design  use.  The  theory  employs 
a  beam-lumped  mass  structural  representation  combined  with  two-dimensional  strip  loading 
modified  for  three-dimensional  flow.21  Figure  9  shows  the  hydroelastic  mode  characteristics 
of  the  AGEH  system  at  29-percent  strut  submergence  as  calculated  by  this  theory.  The  cal¬ 
culation  used  the  only  stiffness  data23  and  the  most  accurate  inertial  data24  available.  A 
slight  simplification  was  introduced  by  using  on  each  subdivision  of  the  submerged  portion 
of  the  structure  constant  steady  loading  va'iies  consisting  of  average  three-dimensional  values. 

Flutter  is  predicted  to  occur  in  tfte  seccnd  hydroelastic  mode  at  a  speed  which  is 
strongly  dependent  on  the  amount  of  structural  damping.  When  the  measured  structural 
damping  of  the  modci  is  added  to  the  calculated  hydrodynamic  damping,  a  flutter  speed  of 
1 1 1  knots  is  obtained.  All  predictions  above  50  knots,  of  course,  are  valid  only  in  terms  of 
scaled  model  speeds.  The  prediction,  corresponding  to  the  broken  line  of  Figure  9,  is 
19-percent  nonconservative  with  respect  to  the  scaled  model  result  of  93  knots.  The  solid 
lines  indicate  predictions  for  the  prototype  system  based  on  measured  structural  damping  of 
the  prototype.  The  slightly  higher  prototype  damping  in  Mode  2  raises  the  predicted  flutter 
speed  to  123  knots.  This  prediction  cannot  be  compared  with  experiment  because  the 
slope  of  the  damping  curve  is  not  known.  The  excellent  agreement  between  predicted  and 
measured  damping  at  intermediate  speeds  indicate-  that  the  flutter  mode  is  predicted 
accurately.  This  agreement  does  not  imply  that  the  theory  can  accurately  predict  damping, 
however,  because  the  damping  of  other,  less  complex  systems  has  not  been  accurately  pre¬ 
dicted  at  all  speeds.1,12 

The  above  comparison  contains  two  known  sources  of  inaccuracy: 

1.  The  model  torsional  stiffness  was  lower  than  the  prototype  stiffness,  as  previously 
described,  so  that  the  prediction  was  apparently  made  for  a  more  rigid  system  than  the 
model.23 


20 

Yates,  E.C.,  Jr.,  “Flutter  Prediction  at  Low  Mass-Density  Ratios  with  Application  to  the  Finite-Span  Noncavitating 
Hydrofoil,”  American  Institute  of  Aeronautics  and  Astronautics  Third  Marine  Systems  and  ASW  Meeting  (Apt-May  1968). 

2 'Liu,  Y.-N.  and  P.K.  Besch,  “Hydrofoil  Flutter  Analysis,  Using  a  Modified  Strip  Theory,”  NSRDC  Report  3624  (Jul 
1971). 

yy 

Rowe,  W.S.  and  T.G.B.  Marvin,  "A  Program  of  Theoretical  Research  on  Hydroelastic  Stability,"  The  Boeing  Company, 
Contract  N00014-67-C-0248  (Nov  1968). 

23Clarlt,  D.J.,  “Instrumentation  and  Static  Calibration  of  the  Foil  Systems  of  the  Navy  Hydrofoil  Boat  PLAJNV1EW 
(AG(EHM),”  David  Taylor  Model  Basin  Report  2149  (Jan  1966). 

24“AGEH-1  Weight  Control  Report-Final,”  Lockheed  Shipbuilding  and  Construction  Company  (Jul  1969). 
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2.  The  straight  beam  structural  representation  was  too  simple  for  the  strut  system.  Frequency 
predictions  obtained  by  using  a  beam  analogy  have  invariably  given  torsional  frequencies  that 
were  too  low  for  this  system.17,25 

Attempts  to  develop  a  constraint  condition  or  parameter  distribution  which  would  improve 
the  prediction  without  distorting  the  mode  shape  have  been  unsuccessful.  However,  a  finite 
element  analogy  which  represented  both  the  forward  and  aft  hinge  points  produced  the 
desired  improvement.*  Therefore,  struts  similarly  supported  at  multiple  root  points  should 
be  analyzed  with  finite  element  techniques.  Flutter  calculations  have  not  been  made  with 
this  method. 

The  importance  of  the  above  sources  of  inaccuracy  could  be  determined  by  further 
analysis  and  possibly  improve  the  prediction.  However,  the  fact  remains  that  only  one  data 
point  is  available  for  evaluation  of  the  theory.  It  is  concluded  that  the  accuracy  of  flutter 
predictions  for  T-foils  cannot  be  determined  for  the  general  case  at  present. 

The  above  calculation  illustrates  the  importance  of  structural  damping  measurements. 

If  the  prototype  damping  were  2  instead  of  9  percent,  the  predicted  flutter  speed  would  be 
83  knots,  which  is  4G  knots  or  33  percent  lower  than  th'  present  prediction.  It  is  clear  that 
damping  must  be  accurately  known  if  accurate  predictions  of  torsional  flutter  are  to  be 
eventually  obtained. 


Flutter  Predictions  for  Struts  and  Struts  with 
Pods 


Flutter  speed  predictions  for  torsion-type  struts  and  struts  with  pods  have  a  fairly  well- 
defined  relationship  to  experiment,  based  on  a  substantial  number  of  comparisons.  Although 
inaccurate  in  most  cases,  the  magnitude  and  direction  of  inaccuracy  correlate  well  with  the 
mass  ratio  of  the  system.  This  correlation  enables  flutter  speed  to  be  estimated  from  the 
calculated  value  with  greatly  improved  accuracy.  A  similar  correction  has  not  yet  been 
determined  for  bending-lype  struts. 

All  available  flutter  speed  predictions  for  struts  without  foils  that  employ  the  modified 
strip  theory  deviate  from  experiment  in  the  manner  shown  in  Figure  1 1.  Struts  with  blunt- 
based  profiles  (included  in  Figure  11)  do  not  exceed  5-percent  thickness  and  are  considered 
comparable  to  subcavitating  struts.  Predictions  for  bending- type  struts  ranged  from  17  to 


25 

Ransleben,  G.E.,  Jr.  and  P.A.  Cox,  Jr.,  “Comparison  of  Measured  Natural  Frequencies  of  the  One-Eighth  Scale  AG(EH) 
Flutter  Model  with  Calculated  Full-Scale  Frequencies,"  Southwest  Research  Institute  Project  02-2260,  Contract  N00600-68- 
D-0375  (Mar  1968). 

* 

Private  communication  from  R.E.  Schauer,  NSRDC  Code  1966. 
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Accuracy  of  Modified  Strip  Theory  Flutter  Speed  Predictions 
for  Strut  Systems  without  Foils 


1 24  percent  nonconservative,  except  for  one  result  which  was  erroneously  predicted  as 
torsional  flutter.  These  bending  flutter  boundaries  correspond  to  a  “new  mode”  instability 
and  not  to  the  instability  which  is  erroneously  predicted  to  occur  in  Mode  3.1  Predictions 
for  torsion-type  struts  varied  from  59  percent  conservative  to  36  percent  nonconservative, 
crossing  over  experimental  values  at  approximately  Attorsion  =  1.9.  A  dramatic  shift  in  the 
direction  of  the  deviation  occurred  when  the  flutter  mode  changed  from  bending  to  torsion. 

The  large  number  of  comparisons  for  torsion-type  struts  makes  it  possible  to  derive  a 
correction  to  calculated  flutter  speeds  for  that  class  of  struts.  Comparisons  are  too  few  to 
permit  this  to  be  done  for  bending-type  struts.  Since  many  additional  data  for  bending-type 
struts  are  available,  the  necessary  calculations  should  be  carried  out  to  determine  whether  a 
similar  correction  is  feasible. 

The  calculated  torsional  flutter  speed  of  a  strut  can  be  corrected  as  follows.  First,  the 
flutter  mode  of  the  strut  must  be  determined  to  ensure  that  torsional  flutter  will  occur.  This 
can  be  done,  as  previously  described,  by  examining  the  in-water  mode  shape  of  the  second 
vibration  mode  of  the  strut.  Second,  the  torsional  mass  ratio  ^torJjon  must  be  calculated  for 
the  strut.  Finally,  the  approximate  percentage  correction  can  be  obtained  from  the  curve  in 
Figure  11.  It  is  important  to  note  that  pod  hydrodynamics,  represented  by  slender  body  theory 
(see  page  417  in  Bisplinghoff  et  al.6)  must  be  included  when  using  the  suggested  correction. 

It  is  not  equivalent  to  use  only  the  pod  added  mass. 

Although  bending  flutter  calculations  cannot  be  corrected  by  using  available  data,  it  may 
be  possible  to  obtain  a  more  accurate  prediction  by  using  two-dimensional  loading.  Both 
the  result  shown  in  Figure  1 1  and  independent  calculations  by  Rowe  and  Marvin22  for  the 
strut  at  Mbending  =  0.101  and  0.262  showed  great  improvement  when  two-dimensional  rather 
than  three-dimensional  loading  was  used.  Further  comparisons  would  be  helpful  in  establishing 
this  point. 

The  discrepancy  between  theory  and  experiment  has  been  shown  to  lie  in  the  hydro- 
dynamic  loading  used.22  The  only  suggested  remedy  has  been  to  modify  the  calculated 
pressure  distribution  to  account  for  possible  cavitation  patterns.22  Such  modifications 
strongly  affected  several  flutter  speed  predictions  in  a  favorable  direction.  However,  a 
systematic  approach  is  not  available  for  determining  the  extent  of  cavitation.  Cavitation  is 
known  to  affect  all  profiles  of  practical  thickness  as  speeds  approach  50  knots.  Comparisons 
of  experimental  and  theoretical  results  at  much  lower  speeds,12  however,  show  that  there  are 
other  sources  of  error  in  addition  to  cavitation.  Improved  accuracy  will  require  improvements 
in  hydrodynamic  loading  theory. 
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DIVERGENCE  OF  SUBCAVITATING  STRUT  SYSTEMS 


KNOWN  PARAMETER  TRENDS 
Vulnerable  Systems 

Both  bending-type  and  torsion-type  struts  can  diverge.  Unswept  bending-type  struts 
probably  do  not  flutter,  and,  therefore,  divergence  is  expected  tc  be  the  only  mode  of  hydro¬ 
elastic  instability  for  such  struts. 

Sweep  Angle 

Divergence  speeds  are  strongly  dependent  on  sweep  angle.  Calculations  summarized  by 
Abramson  et  al.26  indicate  that  divergence  will  rot  usually  occur  above  a  sweep  angle  of 
5  degrees.  Divergence  speed  rapidly  decreases  as  a  strut  is  swept  forward. 

Elastic  Axis  Location 

Divergence  speed  increases  as  the  elastic  axis  is  moved  forward  from  a  position  aft  of 
the  aerodynamic  center,  and  it  becomes  infinite  when  the  elastic  axis  coincides  with  the ' 
aerodynamic  center  of  the  submerged  lifting  surface.  Very  accurate  hydrodynamic  character¬ 
istics  must  be  used  to  calculate  divergence  speed  for  elastic  axis  locations  near  quarter  chord.12 

Submergence 

Increasing  the  wetted  area  lowers  the  divergence  speed  of  a  surface-piercing  strut  or  foil 
because  the  magnitude  of  the  destabilizing  hydrodynamic  moment  is  increased.  Therefore, 
divergence  characteristics  should  be  determined  for  the  condition  of  maximum  wetted  area. 

Size  and  Stiffness 

For  subcavitating  systems  with  equal  values  of  all  nondimensional  parameters,  the 
divergence  speeds  of  a  model  and  a  prototype  are  related  as  follows: 


2®Abramion,  H.N.  et  al.,  “Hydroelaaticity  with  Special  Reference'to  Hydrofoil  Craft,”  NSRDC  Report  2557  (Sep  1967). 
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This  relationship  is  the  same  as  that  for  fiutter  speeds  of  model  and  prototype.  Therefore, 
the  flutter  and  divergence  speeds  bear  the  same  relationship  to  each  other  in  a  scaled  model 
as  in  its  prototype.  If  a  model  encounters  flutter  prior  to  divergence,  so  will  the  prototype. 

If  structural  failure  due  to  approaching  divergence  occurs  prior  to  flutter,  flutter  of  the  proto¬ 
type  will  not  occur  below  the  scaled  speed  reached  by  the  model. 

The  above  expression  is  valid  for  scaling  both  experimental  and  theoretical  divergence 
boundaries  when  the  theories  use  the  same  parameters.  However,  experimental  results  are 
generally  limited  in  accuracy  because  models  never  attain  their  divergence  speed.  Instead, 
they  foil  at  a  lower  speed  which  will  generally  not  simulate  the  failure  of  the  prototype. 


Effect  of  Pods 

Pods  should  lower  the  divergence  speed  of  a  strut.  Theoretically  (according  to  slender 
body  theory),  a  pod  exerts  a  pure,  destabilizing  couple  on  the  tip  of  a  strut  without  con¬ 
sidering  the  attached  strut  or  foils;  see  page  417  in  Bisplinghoff  et  al.6  No  calculations  of 
the  magnitude  of  the  effect  have  been  made. 


Effects  of  Foils 

Foils  should  also  lower  the  divergence  speed  of  a  strut.  An  attached  foil  is  partially 
effective  as  a  reflecting  plate,  thereby  increasing  the  lift  and  moment  on  the  strut.  Similarly, 
the  presence  of  a  yawed  strut  has  been  shown  to  affect  the  loading  on  an  attached  foil.27 
No  divergence  calculations  which  consider  this  interaction  are  available. 


27French,  H.,  “Hydrofoil  Patrol  Craft  PC(H)  Fwd  Hydrofoil  Vibration  Modes,”  The  Boeing  Company,  Aerospace  Division, 
Bureau  of  Ships  Plan  PC<H)  -  1  518  1993317  (Oct  1961). 
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IMPLICATIONS  FOR  PROTOTYPE  SYSTEMS 


Several  design  families  of  T-foils,  based  on  existing  T-foil  structures,  are  free  from 
divergence  or  undesirable  static  deflections  in  subcavitating  flow  within  the  size  range  of 
existing  craft.  Model  experimental  results  support  the  use  of  modified  strip  theory  calcu¬ 
lations  for  predicting  divergence  characteristics. 


AGEH  Main  Design  Family 

Divergence  information  for  the  design  family  based  on  the  AGEH  main  T-foil  is  provided 
by  the  previously  described  flutter  test  of  the  1 /8-scale  AGEH  model.  Since  no  instability 
occurred  prior  to  a  model  speed  of  33  knots,  the  full-scale  system  has  a  divergence  speed 
above  93  knots  under  subcavitating  conditions.  Therefore,  the  full-scale  main  T-foil  system 
will  experience  no  undesirable  deflections  up  to  a  speed  of  50  knots.  All  geometrically 
similar  systems  will  also  be  free  from  static  deflections.  Parameters  for  the  existing  full-scale 
system  are  given  in  Tables  7  and  8. 


Other  Design  Families 

Other  operational  T-foils  have  been  free  from  the  large  deflections  which  precede 
divergence.  Therefore  those  existing  T-foils  define  acceptable  design  families  comprised  of 
geometrically  similar  systems.  System  parameters  for  three  additional  such  systems  are 
summarized  in  Tables  7  and  8.  Possible  structural  differences  in  very  large  struts  limit  the 
applicability  of  these  design  families  to  systems  no  larger  than  the  AGEH  main  T-foil. 


Validity  of  Divergence  Calculations 

Calculations  of  divergence  characterises  using  modified  strip  theory  are  expected  to 
give  accurate  design  information.  Confidence  in  the  theory  is  based  on  two  studies  of  strut 
divergence  behavior  which  gave  excellent  agreement  between  theory  and  experiment.12,28 


28Caapir,  J.R.,  "Divergence  Aiulyiis  of  Swept  Hydrofoils-Computer  Program  (SWDIVRG),”  NSRDC  Report  4245 
(Apr  1974). 
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The  strut  models  included  a  streamlined  profile  model  and  a  thin,  blunt-based  profile  model. 
Both  studies  analyzed  model  divergence  characteristics  from  the  standpoint  of  the  large 
deflections  or  stress  levels  that  precede  divergence.  Excessive  deflections  occurred  at  speeds 
ranging  from  33  to  93  percent  of  accurately  predicted  divergence  speeds.  However,  excessive 
deflections  at  small  fractions  of  the  divergence  speed  occurred  only  when  higher  modes 
appeared  in  the  solution.  When  only  the  fundamental  bending  and  torsion  modes  were  in¬ 
volved,  an  operating  limit  of  60  percent  of  the  divergence  speed  was  adequate.  The  anomalous 
behavior  can  be  detected  by  checking  the  mode  shapes  and  deflection  magnitudes  as 
functions  of  speed. 

These  results  indicate  that  lifting-surface  structures  must  be  restricted  to  speeds  far  be¬ 
low  their  divergence  speeds.  Restrictions  must  be  based  on  avoiding  unacceptable  deflections 
and  stress  levels  for  particular  structural  systems. 


DIVERGENCE  THEORY 

Convenient  divergence  analysis  procedures  have  recently  been  developed  by  Caspar.28 
A  computer  program  is  now  available  which  performs  divergence  calculations  for  any  lifting- 
surface  structure.  A  formula  is  given  for  direct  calculation  of  divergence  speed  for  uniform 
cantilever  struts  or  foils. 

The  computer  program  presented  by  Caspar28  is  based  on  a  strip  loading  theory  that 
has  been  modified  for  three-dimensional  flow  and  combined  with  a  lumped  parameter 
structural  representation.  The  loading  must  be  supplied  separately  by  the  user  from  a  lifting 
surface  theory  or  obtained  from  curves  provided.  The  structural  analogy  for  other  than 
uniform  cantilever  struts  must  be  obtained  from  the  SADSAM  V  structural  analysis  program29 
or  its  equivalent.  In  addition  to  deflections,  the  program  outputs  include  bending  and 
torsional  moments  which  can  be  used  to  determine  stress  levels  for  cantilever  beam  constraint 
conditions. 


MODEL  EXPERIMENTS 

It  does  not  appear  productive  to  use  model  experiments  to  determine  divergence  speeds 
of  structures.  Divergence  theory  will  be  sufficiently  accurate  in  most  cases,  and  models 


29Pe*erson,  L.,  “Theoretical  Basis  for  SADSAM  Computer  Program,”  MacNeal-Schwendler  Corporation  Project  Report 
(Dec  1970). 
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permit  only  approximations  to  be  obtained  because  of  their  predivergence  failure.  If 
experiments  are  performed,  however,  provisions  should  be  made  for  producing  small  static 
deflections  of  the  structure  so  that  a  flow-induced  deflection  boundary  can  be  obtained.  An 
excitation  technique  for  this  purpose  is  described  by  Besch  and  Liu.12  Further  experimental 
experience  in  divergence  testing  is  described  by  Hil borne'  and  by  Squires.9  The  “load 
boundaries”  encountered  by  Squires  are  correlated  with  divergence  speeds  by  Caspar.28 


FLUTTER  OF  CAVITATING  STRUT  SYSTEMS 

Cavitation  and  ventilation  lower  the  flutter  speed  of  at  least  some  lifting  surfaces.  When 
the  higher  speeds  associated  with  cavitation  are  taken  into  account,  cavitating  and  ventilated 
systems  are  potentially  much  more  susceptible  to  flutter  instability  than  are  subcavitating 
systems. 

This  section  discusses  the  fl  fferences  between  the  flutter  of  cavitating  and  /entilated 
systems  and  that  of  subcavitating  systems.  These  differences  consist  of  changes  in  flutter 
speed  produced  by  ventilation  and  the  need  for  an  additional  nondimensional  parameter 
when  flutter  boundaries  are  scaled  in  the  presence  of  cavitation  and  ventilation. 

Also  discussed  are  the  design  implications  of  T-foil  models  and  a  prototype  system. 
Design  guidance  is  available  from  the  prototype  system  up  to  a  speed  of  62  knots. 


THREE  TYPES  OF  CAVITIES 

Cavitation  is  a  local  flow  condition  which  affects  the  hydrodynamic  loading  on  a  sur¬ 
face.  Three  types  of  cavities  occur  in  operating  hydrofoil  systems:  vapor  cavitation,  natural 
ventilation,  and  artificial  ventilation.  Vapor  cavitation  refers  to  submerged  cavities  in  which 
the  pressure  is  near  that  of  water  vapor.  Natural  ventilation  occurs  when  cavities  are  filled 
with  air  at  atmospheric  pressure,  usually  because  the  cavity  is  directly  connected  to  the 
atmosphere.  Artificial  ventilation  corresponds  to  cavities  into  which  air  is  injected. 

Parameters  which  determine  the  occurrence  and  distribution  of  cavitation  and  ventilation 
are  the  geometry  of  the  surface  in  contact  with  the  flow  and  flow  conditions.  Nondi- 
mensionally,  the  system  shape,  cavitation  number  oc,  and  Froude  number  Fn  describe 
system  cavitation  characteristics.  Comparatively  speaking,  vapor  cavitation  corresponds  to 
high  values  of  oc,  artifical  ventilation  to  intermediate  values,  and  natural  ventilation  to  low 
values.  These  variations  in  oc  lead  to  different  scaling  laws  for  each  of  the  three  types  of 
cavitation. 
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KNOWN  PARAMETER  TRENDS 


Lifting-Surface  Profile 

Certain  cavitation  patterns  have  destabilizing  effects  on  lifting  surfaces,  and  profiles  which 
produce  these  cavitation  patterns  should  be  viewed  with  special  caution.  Specifically,  lifting 
surfaces  can  be  greatly  destabilized  by  long  cavities  which  separate  from  the  leading  edge  on 
one  or  both  sides  of  the  surface.  Stability  increases  as  the  separation  point  of  a  long  cavity 
on  one  side  of  a  section  moves  aft  from  the  leading  edge  to  10  percent  of  the  chord. 

At  moderate  speeds,  these  cavitation  patterns  would  be  produced  by  square-nosed  pro¬ 
files,  by  profiles  with  a  spoiler,  and  by  sharp-nosed  profiles  at  sufficiently  large  angles  of 
attack.  At  high  speeds  (above  50  knots),  all  profiles  will  cavitate  aft  of  the  leading  edge 
even  at  small  angles  of  attack,  although  cavity  lengths  might  be  less  than  1  chord.  Since 
nonzero  incidence  angles  are  unavoidable  due  to  sea  conditions,  maneuvering,  or  structural 
oscillation,  lifting  surfaces  operating  at  high  speeds  will  experience  flow  conditions  which  are 
similar  to  those  known  to  be  destabilizing.  For  this  reason,  all  reliable  simulations  of  high¬ 
speed  systems  must  model  the  cavitating  flow  conditions;  stability  estimates  based  on  fully 
wetted  systems  might  be  unconservative. 

The  destabilizing  effect  of  cavitation  was  observed  on  a  surface-piercing  strut1  and  on  a 
two-dimensional  hydrofoil.30,31  The  surface-piercing  strut  data  shown  in  Figures  7  and  8 
exhibited  a  decrease  of  35  percent  in  flutter  speed  when  a  blunt-based  profile  was  changed 
to  a  blunt-nosed  profile  at  a  model  speed  high  enough  to  produce  a  long  cavity.  A  two- 
dimensional  wedge-shaped  hydrofoil  had  a  flutter  speed  as  much  as  34  percent  lower  when  a 
long  rather  than  a  short  (1/2-chord)  cavity  separated  from  the  leading  edge.30  There  was 
little  difference  in  flutter  speed  between  a  wedge  foil  with  a  short,  vapor  cavity30  and  a 
stieamlined  foil.31  This  suggests  that  only  long  cavities  have  a  destabilizing  effect,  at  least 
for  some  foils.  Increases  in  stability  as  the  cavity  separation  point  moves  aft  from  the 
leading  edge  have  been  both  calculated  and  observed  for  two-dimensional  hydrofoils.32 


30Besch,  P.K.,  “Flutter  and  Cavity-Induced  Oscillation  of  a  Two-Degiee-of-Freedom  Hydrofoil  in  Two-Dimensional 
Cavitating  Flow,"  NSRDC  Report  3000  (Apr  1969). 

31Cieslowski  D.S.  and  P.K.  Besch,  “Flutter  of  a  Two-Degree-of-Freedom  Hydrofoil  in  Two-Dimensional  Subcavitating 
Flow,”  NSRDC  Report  3183  (Jan  1970). 

32Song,  C.S.,  “Flutter  of  Supercavitating  Hydrofoils-Comparison  of  Theory  and  Experiment,”  Journal  of  Ship  Research, 
Vol.  16,  No.  3,  pp.  153-166  (Sep  1972). 
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Size  and  Stiffness 


The  characteristics  of  a  cavitating  or  ventilated  flutter  model  have  a  known  relationship 
to  flutter  characteristics  of  prototype  systems  of  different  sizes.  Unlike  subcavitating  proto¬ 
types,  however,  cavitating  prototypes  generally  cannot  be  of  arbitrary  stiffness.  In  some 
cases  at  least,  the  stiffnesses  El  and  GJ  of  the  prototype  must  bear  specific  relationships  to 
the  stiffnesses  of  the  model.  This  restriction  is  created  by  the  need  to  maintain  equal 
values  of  cavitatior  number  oc  at  equivalent  depths  of  model  and  prototype.  Velocity  and 
size  scaling  laws  are  given  below  for  each  type  of  cavitation,  along  with  related  stiffness 
requirements. 

Prototype  and  model  flutter  speeds  and  stiffnesses  must  satisfy  the  following  general 
relationships  for  all  types  of  cavitation: 


U 


f,  prototype 


U 


f,  model 


C  model  / 

5  V1 

prototype  J 


El  O 

prototype  K  model 


El  o 

model  *  prototype 


or 


model 


i 


^  prototype  P  model 


^prototype  V  GJmodel  prototype 


with  the  requirement  that 


El 


C  prototyped  +  P&  "  Pc> 


c 'prototype 


prototype 


C  model  +  P^1  Pc^model 


El 


model 


GJ 


prototype 


prototype  ^Pa  +  P®^1  Pc  ^  prototype 

- -  model 

c  model  ^ Pa  P®^1  —  Pc^model 


For  some  types  of  cavitation,  these  relationships  can  be  simplified  in  the  following  manner. 


Vapor  Cavitation.  Vapor  cavities  have  no  air  inflow;  cavity  pressures  pc  are  very  low,  only 
slightly  above  the  pressure  of  water  vapor  (0.4  psi).  When  both  model  and  prototype 
operate  beneath  an  atmospheric  pressure  of  14.7  psi  at  depths  much  less  than  34  ft 
(corresponding  to  a  pressure  head  much  less  than  atmospheric  pressure),  the  stiffness 
relationships  become 


El  =s 

prototype 


prototype 


model 


El 


model 
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GJ 


prototype 


prototype  . 

C4 


GJ 


model 


model 


Under  these  conditions,  the  flutter  speeds  of  model  and  prototype  are  approximately  equal 


f,  prototype 
f,  model 


/^model 
V  ^  prototype 


a  1.0 


Natural  Ventilation.  Naturally  ventilated  cavities  are  defined  as  having  a  cavity  pressure  pc 
equal  to  atmospheric  pressure.  This  will  occur  for  cavities  with  a  wide  opening  to  the  at¬ 
mosphere.  When  pc  =  p  ,  stiffness  relationships  become 


El 


prototype 


prototype  , 
C3 


El 


model 


model 


GJ 


C  prototype 


prototype  , 

C3 


GJ 


model 


model 


The  resulting  flutter  speed  ratio  is 


U 


f,  prototype 


u 


f,  model 


4 


P  model  prototype 


^prototype  CmodeI 


prototype 


model 


High  speeds  and  restricted  air  paths  resi.it  in  cavity  pressures  below  atmospheric  pressure. 
Such  cavity  conditions  must  be  treated  in  the  same  manner  as  artificially  ventilated  cavities. 


Artificial  Ventilation.  Artificially  ventilated  cavities  have  a  wide  range  of  cavity  pressure 
values,  depending  on  the  air  flow  rate  and  other  conditions.  The  general  scaling  relation¬ 
ships  given  above  must  be  used  to  account  for  these  wide  pressure  variations. 
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In  addition  to  the  above  categories,  combinations  of  two  or  three  types  of  cavitation 
can  occur.  Such  cases  would  have  to  be  scaled  approximately  according  to  the  cavitation 
considered  most  important.  This  is  expected  to  be  the  type  of  cavitation  that  covers  the 
largest  area  on  the  system  lifting  surfaces. 

Under  certain  cavitation  conditions,  it  is  probable  that  the  above  scaling  laws  need  not 
be  fulfilled  and  that  subcavitating  scaling  can  be  used.  An  example  is  the  existence  of 
cavities  much  longer  than  1  chord;  a  change  in  cavity  length  (cavitation  number)  has  little 
effect  on  hydrodynamic  loading.  When  ihe  loading  is  independent  of  cavitation  number,  the 
latter  can  be  neglected  and  the  above  restrictions  on  stiffness  are  not  needed. 


Effects  of  Foils 

A  thin,  blunt-based  profile  strut-pod  system  can  either  be  stabilized  or  destabilized  by 
the  addition  of  a  relatively  large,  lightweight  foil.  A  foil  is  believed  to  change  the 
generalized  mass  ratio  of  the  strut  system.  The  resulting  flutter  speed  changes  may  be 
qualitatively  similar  to  the  curve  in  Figure  5. 

Huang8  found  that  when  a  foil  was  added,  the  flutter  speed  of  a  naturally  ventilated 
strut-pod  system  could  increase  as  much  as  75  percent  or  decrease  by  as  much  as  18  percent 
as  a  function  of  strut  submergence  C.  These  data  are  not  considered  quantitatively  reliable 
for  scaling  purposes  because  the  strut  was  near  a  buckling  condition,  because  the  strut  had 
a  thickness  of  only  3.6  percent,  these  trends  are  probably  accurate  for  subcavitating  struts 
as  well. 


IMPLICATIONS  FOR  PROTOTYPE  SYSTEMS 

Much  less  design  information  is  available  for  cavitating  T-foils  than  for  fully  wetted  T- 
foils.  Only  one  cavitating  T-foil  design  family  is  available  for  design  guidance,  and  this  is 
limited  to  a  known  stable  speed  of  62  knots.  Other  cavitating  designs  are  stable  to  higher 
speeds  but  have  not  been  operated  as  full-sized  prototypes.  These  results  are  discussed 
below. 


DENISON  Tail  Design  Family 

The  DENISON  tail  T-foil,4  partially  described  in  Tables  12  and  13,  defines  a  design 
family  of  geometrically  scaled  T-foils  which  are  known  to  be  stable  up  to  a  speed  of  62  knots. 
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TABLE  12  -  DIMENSIONED  PARAM^  .RS  FOR  FULL-SCALE  AND  MODEL  T.FOILS 
WITH  BLUNT-BASED  STRUT  PROFILES 


Strut 

FRESH-1 

1/4- Scale 
Transiting 

Tail 

FRESH-1 

1/4-Scale 

Transiting 

Main 

FRFSH-1 

Demonstration 

Outboard 

FRESH-1 

Demonstration 

Centerline 

DENISON 

Tail 

L,  in. 

35.7 

28.8 

113 

113 

261.3 

c  (3/4  span)  in. 

5.84 

7.5 

28.8 

28.8 

94 

A  (1/4  chord),  deg 

9 

9 

8.0 

8.0 

0.0 

xea  W*  span>-  in- 
Xcq  ^/4  spanL  in- 
m  (3/4  span),  Ib/in. 

1  (3/4  span),  lb-in. 

El  (3/4  span),  lb-in.2 

4.4 

2.66 

0.402 

2.40 

2.4  x  106 

6.0 

3.06 

0.491 

9.68 

4.3  x  106 

Not  ai 

/ailable 

Not  av 

ailable 

Not  a\ 

ailable 

GJ  (3/4  span),  lb-in.2 

2.3  x  106 

5.1  x  106 

Not  available 

Not  available 

Not  av 

ailable 

Profile 

Blunt  based 

Blunt  based 

Parabolic 

Parabolic 

Blunt  based 

Thickness  (3/4  span) 

20  percent 

16  percent 

11  percent 

11  percent 

7.5  percent 

Dihedral  Angle,  deg 

0.0 

10 

Not  available 

0.0 

0.0 

Pod 

Length,  in. 

16.6 

22.1 

29.1 

44 

226 

Max  .  Diameter,  in. 

1.3 

2.0 

Not  available 

Not  available 

30 

w in- 

5.5 

10 

0.0 

15 

0.0 

xcg.  in. 

7.9 

11 

Not  available 

Not  available 

Not  available 

M,  lb 

1.2 

2.2 

Not  available 

Not  available 

6,000 

ly.  Ib-in.2 

39.5* 

10* 

Not  available 

Not  available 

Not  available 
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TABLE  12  (Continued) 


Foil 

FRESH-1 

1/4- Scale 
Transiting 

Tail 

FRESH-1 

1/4- Scale 
Transiting 

Main 

FRESH-1 

Demonstration 

Outboard 

FRESH-1 

Demonstration 

Centerline 

— 

DENISON 

Tail 

Area,  ft2 

0.385 

0.864 

7.46 

7.46 

49.4 

Root  Chord,  in. 

9.2 

13.7 

29.1 

29.1 

82 

Tip  Chord,  in. 

0.87 

1.3 

8.73 

8.73 

31 

Full  Span,  in. 

11.1 

16.6 

56.8 

56.8 

126 

A.  deg 

48.4 

48.4 

18.0 

18.0 

11.4 

x,oi),  in. 

4.5 

2.75 

2 

2 

-13 

M,  lb 

3.16 

10.9 

Not  available 

Not  available 

Not  available 

x  ,  in. 
eg' 

8.9 

11.8 

Not  available 

Not  available 

Not  available 

Profile 

NACA  16-X04.1 

x  =  0.171 

a  =  1.0 

NACA  16X04.1 

x  =  0.171 

a  =  1.0 

Parabolic 

cambered 

Parabolic 

cambered 

NACA  16-X04 

x  =  0.55 

a  =  1.0 

Flutter  Condition 

C,  in. 

3.9 

11.6 

Not  available 

Not  available 

57.5  to  146.5 

Uf,  knots 

89*  * 

94*** 

above  80 

above  80 

above  62 

ff.  Hz 

17.0** 

7.5*** 

Not  available 

Not  available 

Not  available 

Mode  Shape 

Uncertain 

Uncertain 

Not  available 

Not  available 

Not  available 

- 

0.10 

0.06 

Not  available 

Not  available 

Not  available 

Divergence  Condition 

C,  in. 

6.6 

18.9 

Not  available 

Not  available 

Not  available 

UD,  knots 

98.5 

80 

above  80 

above  80 

above  62 

About  vertical  axis  through  strut  centar  of  gravity. 

•  • 

Minimum  flutter  speed  encountered:  foil  angle  of  attack  was  3  deg. 

*  *  * 

Minimum  flutter  speed  encountered;  foil  angle  of  attack  mas  2  deg. 


TABLE  13  -  NONDIMENSIONAL  PARAMETERS  FOR  FULL-SCALE  AND  MODEL 
T-FOILS  WITH  BLUNT-BASED  STRUT  PROFILES 


Strut 

FRESH-1 

1/4-Scale 

Transiting 

Tail 

FRESH-1 

1/4-Scale 

Transiting 

Main 

FRESH-1 

Demonstrat:on 

Outboard 

FRESH-1 

Demonstration 

Centerline 

DENISON 

Tail 

L/c 

6.11 

3.84 

3.92 

3.92 

2.79 

AR  (to  keel) 

3.10 

2.40 

2.77 

2.62 

1.52 

K 

0.026 

0.041 

0.036 

0.036 

0.0 

xga  (3/4  span)/c 

0.80 

0.80 

Not  available 

Not  availa 

Not  available 

Xcg  <3/4  span)/c 

0.49 

0.41 

Not  available 

Not  available 

Not  available 

4m  (3/4  span) /it p  c2 

0.479 

0.308 

Not  available 

Not  available 

Not  available 

16  lmy  (3/4  span)/trp  c4 

0.291 

0.436 

Not  available 

Not  available 

Not  available 

Pod 

Fineness  Ratio 

12.8:1 

11.1:1 

Not  available 

Not  available 

7.5:1 

Diameter/c 

0.223 

0.267 

Not  available 

Not  available 

0.319 

Xno*/c 

0.942 

1.33 

0.0 

0.521 

0.0 

Xcg,p<x/C 

1.35 

1.47 

Not  available 

Not  available 

Not  available 

Mpod/Msmn 

0.065 

0.119 

Not  available 

Not  available 

above  0.52 

'pod^strut 

0.164 

0.360 

Not  available 

Not  available 

Not  available 

Foil 

Croot/c 

1.69 

1.83 

1.0 

1.0 

0.87 

AR 

2.20 

2.20 

3.0 

3.0 

223 

A,  deg 

48.4 

48.4 

18.0 

18.0 

11.4 

T 

0.0952 

0.0950 

0.30 

'  0.30 

0.38 

Xfoil/c 

0.77 

0.37 

0.069 

0.069 

-0.14 

Xcg/C 

1.52 

1.57 

Not  available 

Not  available 

Not  available 

Mfoi,/MS„u« 

0.168 

0.602 

Not  available 

Not  available 

Not  available 
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TABLE  13  (Continued) 


Flutt-  Condition 

FRESH-1 

1 /4-Scale 
Transiting 
Tail 

FRESH-1 

1/4-Scale 

Transiting 

Main 

FRESH-1 

Demonstration 

Outboard 

FRESH-1 

Demonstration 

Centerline 

DENISON 

Tail 

AR  (strut  wetted  are:! 

0.69 

1.51 

0.80 

0.73 

Not  available 

knf 

0.17 

0.093 

Not  available 

Nov  available 

Not  available 

4  El  (3/4  span)/trp  c2  ;o,2  L4 

4.2  x  10-3 

5.6  x  10~2 

Not  available 

Not  available 

Not  available 

16  GJ  (3/4  span)/irp  L2  u>f2  c4 

0.73 

3.9 

Not  available 

Not  available 

Not  available 

Rnf  (3/4  span) 

8.0  x  106 

8.1  x  106 

above  2.7  x  1C7 

above  2.7  x  107 

above  5.6  x  107 

Fn, 

38 

35 

above  15 

above  15 

above  5.3 

af  (at  surface) 

0.092 

0.083 

below  0.11 

below  C.1 1 

below  0.19 

Divergence  Condition 

AR  (strut  wetted  area) 

1.16 

2.40 

0.80 

0.73 

Not  available 

El  (3/4  span)/2  trp  UD2  L4 

6.4  x  10~4 

4.1  x  10-3 

Not  available 

Not  available 

Not  available 

4  GJ  (3/4  span)/2  ffp  UD2  c2  L2 

0.11 

0.29 

Not  available 

Not  available 

Not  available 

Rn0  (3/4  span) 

6.7  x  106 

6.9  x  106 

above  2.7  x  107 

above  2.7  x  107 

above  5.6  x  107 

FnD 

42 

30 

above  15 

above  15 

above  5.3 

aD  (at  surface) 

0.075 

0.11 

below  0.11 

below  0.1 1 

below  0.19 

59 


As  previously  described,  all  members  of  a  design  family  have  components  (skins,  ribs,  spars, 
etc)  that  are  geometrically  scaled  from  the  parent  and  material  properties  equal  to  those  of 
the  parent.  This  family  should  not  be  extrapolated  to  larger  systems  without  considering 
possible  structural  differences. 


FRESH -I  Transiting  T-Foils 

Quarter-scale  flutter  models  of  the  FRESH-I  transiting  main  and  tail  T-foils  were  tested 
by  Mitchell  and  Rauch2  in  order  to  determine  the  flutter  speeds  of  the  prototype  systems. 

The  transiting  T-foils  were  constructed  in  1963  as  a  second-generation  set  of  foils  for 
FRESH-I,  but  these  have  never  been  used  on  the  craft.  Model  parameters  are  listed  in 
Tables  12  and  13.  Judging  by  the  low  reduced  frequencies  of  flutter,  both  T-foils  experienced 
bending  flutter. 

The  model  flutter  results  can  be  scaled  to  the  size  of  both  the  FRESH-1  prototype  and 
larger  prototypes.  Cavitation  on  the  models  consisted  of  vapor  cavitation  0"  the  foils  and 
natural  ventilation  aft  of  the  blunt-based  struts,  although  cavity  pressure  may  have  been 
below  atmospheric  at  the  high  speed  involved.  Although  cavitation  number  a  would  not  be 
preserved  for  systems  of  different  sizes,  the  cavities  were  long  and  therefore  loading  would  be 
insensitive  to  differences  in  a.  No  observations  were  reported  of  cavitation  separating  from 
the  leading  edge  of  the  struts.  Since  cavitation  number  can  be  neglected,  it  is  possible  to 
scale  flutter  speeds  to  systems  of  different  size  and  arbitrary  stiffness  by  using  subcavitating 
scaling  relationships. 

The  full-scale  transiting  systems  are  geometrically  identical  in  all  parts  an  '  material  to 
the  models,  and  therefore  their  stiffnesses  are  higher  by  the  fourth  power  of  the  length  ratio. 
This  relationship  yields  equal  flutter  speeds  for  model  and  full-scale  systems  since  the  models 
were  tested  in  seawater.  The  full-scale  main  and  tail  systems  are  therefore  expected  to 
flutter  at  94  and  89  knots,  respectively. 

Scaling  these  model  results  to  larger  prototypes  than  the  FRESH-I  systems  requires 
knowledge  of  prototype  stiffnesses.  It  is  unlikely  that  large  prototypes  geometrically  similar 
to  the  transiting  systems  would  be  built  because  such  prototypes  would  be  more  rigid  than 
any  of  the  other  systems  shown  in  Figure  10.  The  PCH-PHM  design  is  the  most  rigid  design 
yet  produced  for  craft  above  100  tons,  and  it  probably  represents  a  practical  maximum  for 
stiffness. 

The  bending  stiffness  of  the  PCH  strut  gives  flutter  speeds  of  82  and  63  knots, 
for  T-foils  similar  to  the  FRESH-I  transiting  main  and  tail  systems,  respectively.  This  result 
applies  to  systems  on  craft  of  100-  to  200-ton  displacements.  It  is  not  known  whether 
larger  designs  could  maintain  this  degree  of  stiffness. 


These  projections  along  with  the  divergence  speed  projections  described  in  another 
section  indicate  that  the  FRESH-1  transiting  T-foil  designs  are  not  acceptable  for  full-sized 
operational  craft  at  speeds  near  80  knots.  It  is  not  clear  which  parameters  are  responsible 
for  this  degree  of  instability.  Structurally,  the  transiting  systems  differ  from  subcavitating 
T-foils  most  markedly  in  their  tapered  struts  and,  according  to  the  description  given  by 
Squires,33  in  the  far  aft  location  of  the  elastic  axis  at  80  percent  of  the  chord. 


FRESH-I  Demonstration  T-Foils 

Two  additional  T-foil  designs,  the  FRESH-1  demonstration  outboard  and  centerline  T- 
foils,  have  been  operated  successfully  under  cavitating  conditions.  However,  these  T-foils 
may  be  of  only  limited  use  for  design  guidance.  An  incomplete  list  of  parameters  for  these 
foils  is  included  in  Tables  12  and  13.  The  FRESH-1  demonstration  foils  are  stable  to  at 
least  80  knots,  but  they  are  small  in  size.  If  these  foil  systems  are  extraordinarily  rigid,  in 
the  manner  of  the  transiting  foils,  and  geometric  scaling  to  larger  sizes  is  impractical, 
reduced  stiffnesses  may  lower  the  known  stable  speed  to  below  80  knots. 


EXPERIMENTAL  MODELING  TECHNIQUES 

Strigent  demands  are  placed  on  model  design  and  experimentation  by  the  need  to  pre¬ 
serve  cavitation  number.  Fortunately,  cavitation  number  must  be  preserved  for  accurate 
simulation  only  when  cavitation  or  ventilation  covers  part  of  a  foil  surface  or  extends  only  a 
short  distance  beyond  the  surface.  In  such  cases,  the  model  must  be  tested  in  a  variable- 
pressure  facility.  Very  high  speeds  at  reduced  pressures  can  be  produced  only  in  the  36-inch 
water  tunnel  at  NSRDC,18  where  a  free  surface  cannot  be  properly  modeled.  Model  design 
must  yield  an  exact  profile  and  stiffness;  in  most  cases,  this  cannot  be  accomplished  with  a 
solid  cross  section  model  but  requires  a  complex  model  design. 

When  hydrodynamic  loading  is  largely  independent  of  cavitation  number,  only  the 
pattern  of  cavitation  on  the  prototype  has  to  be  produced  on  the  model.  This  permits 
modeling  to  be  carried  out  by  using  the  same  approaches  as  for  subcavitating  systems,  except 
that  sufficient  speed  must  be  attained  to  produce  the  required  cavitation. 


33Squircs,  C.F ,  Jr.,  “Preliminary  Flutter  Analysis  of  the  Design  M-51  Hydrofoil  Craft  (FRESH  1),”  Grumman  Aircraft 
Engineering  Co.. .oration  Report  DA  M51-478.1  (Mar  1963). 
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CAVITATING  FLUTTER  THEORY 


Flatter  theory  for  cavitating  systems  in  three-dimensional  flow  is  rudimentary  at  present. 
Development  has  been  held  back  not  only  by  the  limited  amount  of  data  obtained  from 
realistic  T-foil  systems -the  FRESH-I  transiting  systems2  are  the  only  such  systems  for  which 
data  are  available-but  also  by  the  lack  of  unsteady  hydrodynamic  theory  for  all  but  two 
cavitation  patterns.  These  cavitation  patterns  are  short  cavities  separating  from  the  leading 
edge  of  a  foil34  and  long  cavities  separating  from  a  foil  at  an  arbitrary  distance  aft  of  the 
leading  edge.35  No  flutter  analyses  have  been  performed  for  realistic  structures  subject  to 
these  cavitation  patterns.  No  theory  is  available  for  the  many  other  cavitation  patterns  ex¬ 
pected  on  practical  systems,  such  as  blunt-based  profiles  at  high  speeds. 

There  are  many  indications  that  subcavitating  load  theory  is  not  suitable  for  predicting 
flutter  speeds  of  cavitating  strut  systems  other  than  torsion-type  struts  with  thin,  blunt-based 
profiles.  Rowe  and  Marvin2 2 demonstrated  that  calculated  flutter  speeds  are  extremely 
sensitive  to  pressure  distribution  changes  which  simulate  cavitation  (and  to  no  other 
parameter).  Flutter  predictions  for  a  strut  which  probably  cavitated  at  several  locations  along 
its  chord  were  much  more  nonconservative  than  predictions  for  other  cavitating  profiles,1 ,21 ,22 
as  shown  at  ^bcndlng  =  0.262  in  Figure  1 1.  Squires  obtained  “poor”  correla'ton  between 
theory  and  experiment  when  using  subcavitating  load  theory  to  predict  flutter  character¬ 
istics  of  cavitating  T-foils.33  The  lift  slope  of  a  24-percent-thick,  blunt-based  profile  strut 

% 

was  found  to  be  substantially  higher  than  that  of  a  streamlined  strut  of  12-percent  thickness, 
suggesting  that  different  modifications  to  a  strip  theory  would  be  required  for  a  moderately 
thick  blunt-based  strut  and  a  streamlined  strut. 

The  semiempirical  approach  provides  an  alternative  to  pure  theory.  This  approach 
would  involve  measurement  of  the  unsteady  lift  and  moment  on  a  model  system  in  properly 
simulated  flow.  The  measured  loading  would  be  used  with  a  suitable  structural  representation 
to  calculate  flutter  speeds.  The  accuracy  of  this  method  depends  on  the  accurary  of  the 
measured  loading.  Considerable  expense  would  be  involved  in  obtaining  snanwise  distri¬ 
butions  of  loading.36  Because  the  expense  of  an  accurate  measurement  of  load  distribution 


34 

Steinberg,  11.  and  S.  Karp,  “Unsteady  Mow  Past  Partially  Cavitated  Hydrofoils,"  Proceedings  f  ourth  Naval  Hydro¬ 
dynamics  Symposium  (Aug  1962). 

35Kaplan,  P.  and  C.J.  Henry,  “A  Study  of  the  llydroclastic  Instabilities  of  Supercavitating  Hydrofoils,”  Journal  of  Ship 
Research,  Vol.  4,  No.  3,  pp.  28-38  (Dec  I960). 

36Ranslebcn,  (l.f.,  Jr.,  “experimental  Determination  of  Steady  and  Unsteady  Loads  on  a  Surface  Piercing,  Ventilated 
Hydrofoil,"  Journal  of  Ship  Research,  Vol.  13,  No.  I,  pp.  I- II  (Mar  1969). 

* 

Reported  by  N.L.  Dailey  in  an  NSRDC  report  for  internal  use. 


62 


would  probably  be  greater  than  the  expense  of  building  an  elastically  scaled  model  which 
could  be  used  to  establish  flutter  stability  directly,  the  semiempirical  approach  is  of  practical 
use  only  when  it  is  desired  to  avoid  building  a  large  number  of  flutter  models. 

It  is  concluded  that  cavitating  flutter  theory  is  not  now  available  for  general  design  use. 
Based  on  results  in  Figure  1 1,  torsion-type  struts  with  thin,  blunt-based  profiles  can  be 
treated  by  using  subcavitating  flutter  theory. 


DIVERGENCE  OF  CAVITATING  STRUT  SYSTEMS 

Relatively  little  information  on  the  divergence  of  cavitating  systems  is  available  for 
design  use.  Because  of  the  drastic  changes  in  hydrodynamic  loading  which  accompany  cavi¬ 
tation,  cavitating  systems  should  have  divergence  characteristics  which  are  substantially 
different  from  those  of  subcavitating  systems.  Design  is  hindered  by  the  la:k  of  theories  of 
established  accuracy,  leaving  model  experiments  as  the  only  dependable  alternative. 


KNOWN  PARAMETER  TRENDS 
Sweep  Angle 

The  effect  of  sweep  angle  in  raising  divergence  speed  has  been  experimentally  demon¬ 
strated  for  thin,  blunt-based  struts  by  Hilbome7  and  by  Squires.9  Struts  at  less  than 
5  degrees  of  sweep  experienced  large  deflections  and  occasional  yielding,  which  were 
indicative  of  approaching  divergence,  while  deflections  were  insignificant  as  higher  sweep 
angles.  Flutter  replaced  divergence  as  the  mode  of  hydroelastic  instability  at  sweep  angles 
of  5  degrees  and  higher. 


Size  and  Stiffness 

Divergence  speeds  of  cavitating  systems  scale  the  same  as  flutter  speeds  of  cavitating 
systems.  In  some  cases,  the  requirement  to  preserve  cavitation  number  places  restrictions  on 
model  stiffness. 


63 


IMPLICATIONS  FOR  PROTOTYPE  SYSTEMS 


The  DENISON  tail  T-foil,5  described  in  Tables  12  and  13,  establishes  a  design  family  of 
geometrically  scaled,  equivalent  material  systems  which  are  free  from  divergence-related 
deflections  up  to  a  speed  of  62  knots.  Extrapolations  should  not  be  made  to  larger  size 
systems  without  further  study. 

Two  additional  cavitating  T-foil  designs  have  been  found  to  be  subject  to  large  predi¬ 
vergence  deflections,  in  one  case  at  a  speed  below  flutter  inception.  The  designs  are  the  1/4- 
scale  FRESH-I  transiting  main  and  tail  T-foils,2  parameters  for  which  are  given  in  Tables  12 
and  13.  Since  both  systems  were  capable  of  operation  at  speeds  above  the  initially  unstable 
speeds,  the  flow  presumably  changed  in  some  way  at  the  higher  speeds.  It  is  therefore 
obviously  important  to  establish  stability  for  all  possible  operating  conditions,  not  merely  the 
highest  speed. 

According  to  the  scaling  approach  explained  earlier  for  flutter  speeds,  predivergence 
deflections  will  occur  on  the  full-scale  FRESH-I  transiting  main  and  tail  T-foils  at  80  and 
98.5  knots,  respectively,  the  same  speeds  as  for  the  models.  Larger  scale  prototypes  of  both 
the  main  and  tail  systems,  with  stiffnesses  corresponding  to  the  PCH-PHM  design  family 
shown  in  Figure  10,  would  be  affected  at  70  knots.  These  designs  could  not  be  used  with 
confidence  at  speeds  very  far  into  the  cavitating  flow  region. 

Prototype  systems  should  not  be  expected  to  survive  a  static  instability  in  the  manner 
that  the  transiting  system  models  did.  This  result  is  unique  among  experimental  results;  most 
other  instability-related  deflections  have  caused  model  failure  and  no  other  instances  have 
been  successfully  pursued  to  higher  speeds.  The  least  indication  of  instability  in  a  design 
would  of  course  n^t  be  tolerable.  Although  it  is  assumed  that  changes  in  loading  due  to 
varying  amounts  of  cavitation  led  to  survival  of  the  transiting  models,  other  changes,  e.g., 
structural,  could  have  occurred. 


CAVITATING  DIVERGENCE  THEORY 

Divergence  theory  for  cavitating  systems  is  unavailable  for  design  use  because  flutter 
theory  of  cavitating  systems  is  deficient  except  for  thin,  blunt-based  profile  struts.  The 
semiempirical  approach  mentioned  in  discussing  cavitating  flutter  theory  is  a  possible  approach 
for  divergence  predictions  as  well. 
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MODEL  EXPERIMENTS 


Either  solid-section  or  more  complex  models  could  be  used  to  establish  divergence-free 
speed  ranges  of  prototypes.  Care  must  be  taken  to  simulate  the  prototype  cavitation  pattern 
exactly  in  those  cases  where  loading  will  be  sensitive  to  changes  in  cavity  length.  Cavitation 
produced  by  operational  angles  of  attack  must  be  considered  in  the  model  study. 


CONCLUSIONS 

1.  Within  the  size  range  of  existing  hyti.ofoil  craft,  several  design  families  of  T-foils  are 
free  of  flutter  and  divergence  at  all  subcavitating  speeds  or  up  to  approximately  50  knots. 

2.  Within  the  size  range  of  existing  craft,  one  design  family  of  T-foils  is  stable  in  the 
cavitating  speed  range  up  to  a  speed  of  62  knots.  No  stable  designs  nave  been  established 
above  62  knots. 


RECOMMENDATIONS 

1.  The  engineering  aspects  of  constructing  T-foil  systems  larger  than  the  AGEH  main  T- 
foil  should  be  studied  to  determine  the  impact  of  possible  structural  differences  on  hydro¬ 
elastic  stability  and  other  design  criteria. 

2.  Studies  of  flutter  and  divergence  characteristics  should  be  made  for  both  subcavitating 
and  cavitating  T-foils.  Model  testing  and  theoretical  analysis  should  be  used. 

3.  A  technique  foi  measuring  damping  of  full-scale  T-foils  should  be  developed  and 
should  include  installation  and  trial  on  an  existing  craft.  The  technique  should  then  be  used 
on  all  future  craft. 

4.  A  standard  specification  should  be  developed  for  hydroelastic  stability  on  high-speed 

craft. 
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APPENDIX 

KINEMATIC  SCALING  OF  FLUTTER  AND  DIVERGENCE  OF  T-FOILS 


Principles  of  dimensional  analysis  permit  flutter  data  from  one  T-foil  to  be  used  to  pre¬ 
dict  flutter  characteristics  of  other  T-foils  of  different  size  and  stiffness.  Flutter  boundaries 
of  a  prototype  and  its  model  can  be  compared  if  the  two  systems  undergo  equivalent 
motions,  corresponding  to  zero  damping,  when  all  scaling  relationships  are  satisfied.  This 
condition  is  called  kinematic  similarity  or  scaling.37  It  occurs  when  the  ratios  of  displace¬ 
ment,  velocity,  and  acceleration  between  the  two  systems  are  the  same  at  all  points.38 

Scaled  systems  have  equal  values  or  effectively  equal  values  of  all  nondimensional 
parameters.  The  required  nondimensional  parameters  for  kinematic  scaling  can  be  obtained 
either  by  nondimensionalizing  the  equations  of  motion  of  the  system,  as  has  been  done  for  a 
wing  (see  the  material  starting  on  page  699  in  Bisplinghoff  et  al.6)  or  by  forming  dimension¬ 
less  ratios  of  all  parameters  that  should  appear  in  those  equations. 

This  treatment  will  list  the  relevant  parameters  and  the  dimensionless  ratios  obtained 
from  an  unpublished  nondimensionalization  of  the  equations  of  motion.  Some  parameter 
ratios  are  based  on  structural  properties  at  the  3/4-span  position.  These  properties  are 
representative  of  the  overall  properties  of  tapered  structures  (see  page  533  in  Bisplinghoff 
et  al.6). 


SUBCAVITAT1NG  FLUTTER 

Dimensioned  Flutter  Parameters  for 
Subcavitating  T-Foils 

Strut  Parameters:  L,  c,  L  cos  A,  xea  (y),  xcg  (y),  m  (y),  Imy  (y),  El  (y),  GJ  (y),  2 
Pod  Parameters:  length,  maximum  diameter,  xnose,  xcg,  M,  Iy,  1^ 

Foil  Parameters:  L,  c(z'),  L  cos  A,  xfoil,  xea,  (z'),  xcg  (z'),  m  (z'),  ly,  I]nz,  (z'X  1^, 
El  (z'),  GJ  (z') 

System  Parameters:  p,  cs  (in  air) 


37Langhaar,  H.L.,  “Dimensional  Analysis  and  Theory  of  Models,”  John  Wiley  &  Sons,  Inc.,  New  York  (1951). 
38Albertson,  M.L.  et  al.,  “Fluid  Mechanics  for  Engineers,"  Prentice-Hall,  Inc.,  Englewood  Cliffs,  N.J.  (I960). 


Preceding  page  blank 
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Boundary  Conditions  for  T-Foils:  root  support  condition,  strut-foil  connection 
Flutter  Condition:  Uf,  ff,  node  shape 


Nondimensional  Flutter  Parameter  Ratios  for 
Subcavitating  T-Foils 

Strut  Parameters:  L/c,  AR  (to  root),  k,  xm  (y)/c,  x£g  (y)/c,  4m  (y)/rrpc2, 

16  Imy(y) 

’  *eiN  *cgN 

TTpC4 

/’od  Parameters:  fineness  ratio,  (maximum  diameter)/c,  xnoJe/ c,  xcg/c, 

^pod^strut’  ^y, pod ^y .strut’  ^.pod^y.pod 

Foil  Parameters:  cfoot/c,  AR,  A,  t,  xfoa/c,  xM  (z*)/c,  xcg  (z')/c,  mfoiJ  (3/4  span)/mJtrut 
(3/4  span),  Iy,foil/IyiStnit,  Im2..foil  (3/4  span)/Im,/>strut  (3/4  span), 
V,fou/V,Pod>  EIfoa  <3/4  sPan>/E1strut  (3/4  span),  GJfojl  (3/4  span)/ 
strut  (3/4  span),  mN  (H,  Imrt|  (?'),  EIN  (?'),  GJN  (D 

System  Parameters:  fs,  p  (see  text  for  approximate  forms) 


Nondimensional  Boundary  Conditions 


Flutter  Parameters:  AR  (wetted  area),  knf, 
normalized  mode  shape 


4  El  (3/4  span) 
it  pc2  wf2  L4 


16  GJ  (3/4  span) 
irpL2  u}2  c4 


These  parameters  were  obtained  from  a  modified  strip  theory  that  has  been  used  in  a 
large  number  of  hyd  ofoil  flutter  calculations.1  •* 2,21  Some  reduction  was  possible  in  the 
number  of  nondimensional  parameters  because  several  dimensioned  parameters  (L,  c,  A, 

E,  I,  G,  J)  appear  only  in  combination  with  other  parameters. 

Two  additional  dimensionless  ratios,  Reynolds  number  Rn  and  Froude  number  Fn,  are 
theoretically  operative  in  flutter  scaling.  These  parameters  can  usually  be  neglected  for  the 
following  reasons.  Flutter  of  both  model  and  prototype  usually  occurs  at  relatively  high 
speeds  for  which  Rn  is  well  above  4  x  1 05 ,  the  value  above  which  changes  in  Rn  have  little 
effect  on  oscillatory  loading  on  airfoils  (see  page  710  in  Bisplinghoff  et  al.6)  and  should 
similarly  not  !r.fluence  hydrofoil  loading  characteristics.  It  is  therefore  assumed  that  Rn  can 
-be  neglected  for  hydrofoil  flutter. 
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The  effect  of  Froude  number  on  hydrofoil  loading  has  been  studied  for  two-dimensional 
flow  by  Ashley  et  al.39  They  predict  that  as  a  result  of  interaction  between  the  foil  and  the ' 
free  surface,  there  will  be  a  strong  dependence  on  Fn  when  the  foil  is  close  to  the  surface  and 
foil  speed  and  oscillation  frequency  are  low.  It  is  expected  that  this  interaction  will  be 
weaker  in  the  three-dimensional  case.  Since  hydrofoil  systems  will  rarely  if  ever  flutter  at 
sufficiently  low  speeds  and  frequencies,  Froude  number  is  not  important  as  a  flutter 
parameter  in  the  present  case  (subcavitating  flow).  It  is  shown  below,  however,  that  Froude 
number  is  necessarily  preserved  when  natu.-a'ly  ventilated  systems  are  scaled.  Neither  Rn  nor 
Fn  is  required  for  scaling  subcavitating  flutter,  but  both  parameters  are  included  for 
reference  in  experimental  results  given  in  this  report. 


SUBCAVITATING  DIVERGENCE 

The  parameters  used  to  describe  divergence  are  those  flutter  parameters  which  do  not 
involve  oscillation  or  inertia  of  the  structure.  El  and  GJ  are  nondimensionalized  in  the  form 
that  would  appear  in  the  equations  of  motion  of  the  structure. 


Dimensioned  Divergence  Parameters  for 
Subcavitating  T-Foils 

Strut  Parameters:  L,  c,  L  cos  A,  xea  (y),  El  (y),  GJ  (y),  8 

Pod  Parameters:  length,  maximum  diameter,  xnoJe 

Foil  Parameters:  L,  c  (z'),  L  cos  A,  xfoil,  xea  (z'),  El  (z'),  GJ  (z') 

System  Parameter:  p 

Boundary  Conditions  for  T-Foils:  root  support  condition,  strut-foil  connection 
Divergence  Condition:  UD,  mode  shape 


39Ashley,  H.  et  at.,  “New  Directions  in  Lifting  Surface  Theory,”  AIAA  Journal,  Vol.  3,  No.  1,  pp.  3-16  (Jan  1965) 
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Nondimensional  Divergence  Parameter  Ratios 
for  Subcavitating  T- Foils 


Stmt  Parameters:  L/c,  AR  (to  root),  k,  x(t  (y)/c,  EIN  (17),  GJN  (q) 
Pod  Parameters:  fineness  ratio,  (maximum  diameter)/c,  xnofe/ c 
Foil  Parameters:  cf00,/c,  AR,  A,  r,  xfoil/c,  xM  (z')/c,  EIfofl  (3/4  span)/ 
EI«rut  <3/4  sPan)’  GJfoa  <3/4  sPan)/GIftruj  (3/4  span), 
E!n  (“),  GJn  (JO 

Nondimensional  1  Boundary  Conditions 


Divergence  Parameters: 


AR  (wetted  area), 


El  (3/4  span) 


2  jrp  UD2  L4 
normalized  mode  shape 


2  GJ  (3/4  span) 
itp  UD2  c2  L2 


FLUTTER  OF  CAVITATING  AND 
VENTILATED  SYSTEMS 


Scaling  flutter  boundaries  for  cavitating  and  ventilated  systems  requires  that  hydrodynamic 
forces  due  to  the  cavitating  or  ventilated  flow  be  similar  in  model  and  prototype  systems. 

This  will  occur  for  similar  profiles  when  cavitation  number  oc  is  equal  at  all  equivalent 
depths  in  both  systems  where  there  is  cavitation  or  ventilation.  This  may  be  expressed  as 


pa  +  pgli  -  pc 


\_ 

2 


p  U2 


prototype 


pa  +  pgh  -  pc 


model 


cr,  in  terms  of  velocity, 


E  prototype  _  Pmodel  ^Pa  +  P®^1  Pc  ^prototype 
model  P  prototype  (Pg  +  —  Pc)  model 

The  value  of  pc  varies  according  to  the  type  of  cavitation  present.  For  natural  ventilation, 
the  expression  is  greatly  simplified  because,  by  definition,  pc  =  pa.  In  that  case,  Froude 
number  Fn  is  also  preserved  between  the  two  systems.  In  general,  however,  preservation  of 
Fn  in  addition  to  oc  would  place  impractical  requirements  on  pa  and  pc  and  would  not  be 
-  necessary  so  long  as  Fn  is  reasonably  large  (near  5.0). 


It  is  noted  for  reference  that  few  test  facilities  are  capable  of  simulating  cavitation  con¬ 
ditions  that  occur  for  full-scale  prototype  operation  above  50  knots.  Flutter  models  of  high¬ 
speed  prototypes  would  generally  have  to  be  tested  on  the  1 50-knot  carriage  in  the  High- 
Speed  Hydrodynamics  Facility  at  Langley  Field,  Virginia40  or,  if  a  free  surface  is  not 
required,  in  the  36-inch  variable-pressure  water  tunnel21  at  NSRDC  Carderock. 

The  requirement  that  ac  be  preserved  between  model  and  prototype  fixes  the  values  of 
stiffnesses  FI  and  GJ  that  must  be  used  in  a  model.  These  values  vary  according  to  the  type 
of  cavitation  present.  The  scaling  laws  for  each  type  of  cavitation  are  given  in  the  section  of 
this  report  dealing  with  known  parameter  trends  of  cavitating  systems.  Individual  system 
parameters  are  listed  below. 


Dimensioned  Flutter  Parameters  for  Cavitating 
T-Foils 

These  include  pa,  pc.  and  surface  geometry  for  all  submerged  components,  in  addition 
to  all  dimensioned  flutter  parameters  previously  given  for  subcavitating  T-foils. 


Nondimensional  Flutter  Parameter  Ratios  for 
Cavitating  T-Foils 

These  n  elude  oc  and  nondimensionalized  profiles  for  ail  submerged  components,  in 
addition  to  all  nondimensional  flutter  parameters  previously  given  for  subcavitating  T-foils. 


DIVERGENCE  OF  SUBCAVITATING  T-FOILS 

Divergence  parameters  are  the  same  as  flutter  parameters  except  for  the  removal  of 
quantities  involving  structural  inertia  and  oscillation  frequency.  The  parameters  given  below 
were  obtained  from  unpublished  modified  strip  theory  divergence  equations. 


40 
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Dimensioned  Divergence  Parameters  for 
Subcavitating  T- Foils 


Strut  Parameters:  L,  c,  L  cos  A,  xea  (y),  El  (y),  GJ  (y),  8 

Pod  Parameters:  length,  maximum  diameter,  xnoJe 

Foil  Parameters:  L,  c  (z'),  L  cos  A,  xfoi],  xM  (z'),  El  (z'),  GJ  (z') 

System  Parameter:  p 

Boundary’  Conditions  for  T-Foils:  root  support  condition,  strut-foil  connection 

Divergence  Condition:  UD ,  mode  shape 

Nondimensional  Divergence  Parameters  for 
Subcavitating  T-Foils 


Strut  Parameters:  L/c,  AR  (to  root),  «,  xea  (y)/c,  EIN  (r?),  GJN  (tj) 
Pod  Parameters:  fineness  ratio,  (maximum  diameter)/c,  xnosil c 
Foil  Parameters:  k,  AR,  xfojl/c,  xea  (z')/c,  EIN  (O.  GJN  (?') 
Nondimensional  Boundary  Conditions 

El  (3/4  span) 


Divergence  Condition:  AR  (wetted  area), 


itp 

normalized  mode  shape 


2  itp  L4  IL 2  cos2  A 


2  GJ  (3/4  span) 

"  “  > 

itp  c2  L2  UD  ‘  cos2  A 


DIVERGENCE  OF  CAVITATING  AND 
VENTILATED  T-FOILS 


Dimensioned  Divergence  Parameters  for 
Cavitating  T-Fcils 

These  include  pa,  pc,  surface  geometry  for  all  submerged  components,  in  addition  to  all 
dimensioned  divergence  parameters  previously  given  for  subcavitating  T-foils. 


Nondimensional  Divergence  Parameter  Ratios 
for  Cavitating  T-Foils 

These  include  oc  and  nondimensionalized  profiles  for  all  submerged  components,  in 
addition  to  all  nondimensional  divergence  parameters  previously  given  for  subcavitating  T- 
foils. 


72 


REFERENCES 


1.  Besch,  P.K.  and  Y.-N.  Liu,  “Bending  Flutter  and  Torsional  Flutter  of  Flexible 
Hydrofoil  Struts,”  Ninth  Symposium  on  Naval  Hydrodynamics,  Paris,  France  (Aug  1972); 
also  NSRDC  Report  4012  (Feb  1973). 

2.  Mitchell,  L.  and  F.J.  Rauch,  Jr.,  “Dynamic  Tests  of  the  1 /4-Scale  Models  of  the 
80-Knot  Transiting  Strut-Foil  Systems  for  the  FRESH  I  Hydrofoil  Test  Craft,”  Grumman 
Aircraft  Engineering  Corporation,  Contract  NObs  86826  (Aug  1964). 

3.  Kennard,  E.H.,  “Irrotational  Flow  of  Frictionless  Fluids,  Mostly  of  Invariable 
Density,”  David  Taylor  Model  Basin  Report  2299,  pp.  390-392  (Feb  1967). 

4.  Scanlan,  R.H.  and  R.  Rosenbaum,  “Aircraft  Vibration  and  Flutter,”  Dover 
Publications,  Inc.,  New  York  (1968). 

5.  Schaffer,  A.P.,  “Flutter  Flight  Testing  of  a  High  Speed  Hydrofoil  Craft-The  H.S. 
DENISON,”  Grumman  Aircraft  Engineering  Corporation  Report  FT-HYD-Ia.  1  (Apr  1963). 

6.  Bisplinghoff,  R.T.  et  al.,  “Aeroelasticity,”  Addison-Wesley  Publishing  Company,  Inc., 
Cambridge,  Massachusetts  (1955). 

7.  Hilbome,  D.V.,  “The  Hydroelastic  Stability  of  Struts,”  Admiralty  Research 
Laboratory  (Great  Britain)  Report  ARL/R  l/G/HY/5/3  (1958). 

8.  Huang,  T.T.,  “Experimental  Study  of  a  Low  Modulus  Flutter  Model  for  Strut-Foil- 
Pod  Configurations,”  Hydronautics,  Inc.  Technical  Report  459-2  (Jul  1967). 

9.  Squires,  C.E.,  Jr.,  “Hydrofoil  Flutter,  Small  Sweep  Angle  Investigation-Final  Report,” 
Grumman  Aircraft  Engineering  Corporation  Report  DA  Nonr-3989.3  (Nov  1963). 

10.  Baird,  E.F.  et  al.,  “Investigation  of  Hydrofoil  Flutter-Final  Report,”  Grumman 
Aircraft  Engineering  Corporation  Report  DA  10-480-3  (Feb  1962). 

1 1.  Abramson,  H.N.  and  G.E.  Ransleben,  Jr.,  “An  Experimental  Investigation  of 
Fluttei  of  a  Fully  Submerged  Subcavitating  Hydrofoil,”  J.  Aircraft,  Vol.  2,  No.  5,  pp.  439- 
442  (Scp-Oct  1965). 

12.  Besch,  P.K.  and  Liu,  Y.-N.,  “Flutter  and  Divergence  Characteristics  of  Four  Low 
Mass  Ratio  Hydrofoils,”  NSRDC  Report  3410  (Jan  1971). 

1 3.  Woolston,  D.S.  and  G.E.  Castile,  “Some  Effects  of  Variations  in  Several  Parameters 
Including  Fluid  Density  on  the  Flutter  Speed  of  Light  Uniform  Cantilever  Wings,”  NACA 

-IN  2558  (1951). 


73 


14.  Ransleben,  G.E.,  Jr.,  “Experimental  Determination  of  Variation  of  Hydrofoil 
Flutter  Speed  with  Mass  Ratio,”  Southwest  Research  Institute,  Contract  N00014-69-C-0219 
(Apr  1970). 

15.  Ransleben,  G.E.,  Jr.,  “Description  of  a  One-Eighth  Scale  Flutter  Model  of  the  Main 
Hydrofoil  Assembly  of  the  AG(EH)  800  Hydrofoil  Research  Ship,”  Southwest  Research 
Institute  Project  02-1960  (Dec  1966). 

16.  Peoples,  J.R.,  “Frequencies  and  Damping  of  Full  Scale  Hydrofoils  by  “Pluck  Test” 
Methods,”  43rd  Shock  and  Vibration  Bulletin,  Shock  and  Vibration  Information  Center, 

Naval  Research  Laboratory  (Jun  1973). 

17.  “PCH-1  Mod  1  Detail  Design -Structural  Summary,”  The  Boeing  Company  (1969). 

18.  Brownell,  W.F.  and  M.L.  Miller,  “Hydromechanics  Cavitation  Research  Facilities  and 
Techniques  in  Use  at  the  David  Taylor  Model  Basin,”  Symposium  on  Cavitation  Research 
Facilities  and  Techniques,  American  Society  of  Mechanical  Engineers  (May  1964);  ako  David 
Taylcr  Model  Basin  Report  1856  (Oct  1964). 

19.  White,  R.G.,  “Use  of  Transient  Excitation  in  the  Dynamic  Analysis  of  Structures,” 
The  Aeronautical  Journal  of  the  Royal  Aeronautical  Society,  Vol  73,  pp.  1047-1050  (Dec 
1969). 

20.  Yates,  E.C.,  Jr.,  “Flutter  Prediction  at  Low  Mass-Density  Ratios  with  Application  to 
the  Finite-Span  Noncavitating  Hydrofoil,”  American  Institute  of  Aeronautics  and  Astronautics 
Third  Marine  Systems  and  ASW  Meeting  (Apr-May  1968). 

21.  Liu,  Y.-N.  and  P.K.  Besch,  “Hydrofoil  Flutter  Analysis,  Using  a  Modified  Strip 
Theory,”  NSRDC  Report  3624  (Jul  1971). 

22.  Rowe,  W.S.  and  T.G.B.  Marvin,  “A  Program  of  Theoretical  Research  on  Hydroelastic 
Stability,”  The  Boeing  Company,  Contract  N00014-67-C-0248  (Nov  1968). 

23.  Clark,  D.J.,  “lnstrumentat;on  and  Static  Calibration  of  the  Foil  Systems  of  the  Navy 
Hydrofoil  Boat  PLAINV1EW  (AG(EH)-l),”  David  Taylor  Model  Basin  Report  2149  (Jan 
1966). 

24.  “AGEH-1  Weight  Control  Report-Final,”  Lockheed  Shipbuilding  and  Construction 
Company  (Jul  1969). 


74 


25.  Ransleben,  G.E.,  Jr.  and  P.A.  Cox,  Jr.,  “Comparison  of  Measured  Natural  Frequencies 
of  the  One-Eighth  Scale  AG(EH)  Flutter  Model  with  Calculated  Full-Scale  Frequencies,” 
Southwest  Research  Institute  Project  02-2260,  Contract  N00600-68-D-0375  (Mar  1968). 

26.  Abramson,  H.N.  et  al.,  “Hydroelasticity  with  Special  Reference  to  Hydrofoil  Craft,” 
NSRDC  Report  2557  (Sep  1967). 

27.  French,  H.,  “Hydrofoil  Patrol  Craft  PC(H)  Fwd  Hydrofoil  Vibration  Modes,” 

The  Boeing  Company,  Aerospace  Division,  Bureau  of  Ships  Plan  PC(H)  -  1  518  1993317 
(Oct  1961). 

28.  Caspar,  J.R.,  “Divergence  Analysis  of  Swept  Hydrofoils-Computer  Piogram 
<  (SWDIVRG),”  NSRDC  Report  4245  (Apr  1974). 

29.  Peterson,  L.,  “Theoretical  Basis  for  SADSAM  Computer  Program,”  MacNeal- 
Schwendler  Corporation  Project  Report  (Dec  1970). 

30.  Besch,  P.K.,  “Flutter  and  Cavity-Induced  Oscillation  of  a  Two-Degree-of-Freedom 
Hydrofoil  in  Two  Dimensional  Cavitating  Flow,”  NSRDC  Repoit  3000  (Apr  1969). 

31.  Cieslowski,  D.S.  and  P.K.  Besch,  “Flutter  of  a  Two-Degree-of-Freedom  Hydrofoil  in 
Two-Dimensional  Subcavitating  Flow,”  NSRDC  Report  3183  (Jan  1970). 

32.  Song,  C.S.,  “Flutter  of  Supercavitating  Hydrofoils-Comparison  of  Theory  and 
Experiment,”  Journal  of  Ship  Research,  Vol.  16,  No.  3,  pp.  153-166  (Sep  1972). 

33.  Squires,  C.E.,  Jr.,  “Preliminary  Flutter  Analysis  of  the  Design  M-51  Hydrofoil 
Craft  (FRESH  I),”  Grumman  Aircraft  Engineering  Corporation  Report  DA  M 5 1-478. 1  (Mar 
1963). 

34.  Steinberg,  H.  and  S.  Karp,  “Unsteady  Flow  Past  Partially  Cavitated  Hydrofoils,” 
Proceedings  Fourth  Naval  Hydrodynamics  Symposium  (Aug  1962). 

35.  Kaplan,  P.  and  C.J.  Henry,  “A  Study  of  the  Hydroelastic  Instabilities  of  Super¬ 
cavitating  Hydrofoils,”  Journal  of  Ship  Research,  Vol.  4,  No.  3  pp.  28-38  (Dec  1960). 

36.  Ransleben,  G.E.,  Jr.,  “Experimental  Determination  of  Steady  and  Unsteady  Loads 
on  a  Surface  Piercing,  Ventilated  Hydrofoil,”  Journal  of  Ship  Research,  Vol.  13,  No.  1, 
pp.  1-11  (Mar  1969). 

37.  Langhaar,  H.L.,  “Dimensional  Analysis  and  Theory  of  Models,”  John  Wiley  &  Sons, 
Inc.,  New  York  (1951). 


* 


75 


38.  Albertson,  M.L.  et  al.,  "Fluid  Mechanics  for  Engineers,”  Prentice-Hall,  Inc., 
Englewood  Cliffs,  N.J.  (1960). 

39.  Ashley,  H.  et  al.,  “New  Directions  in  Lifting  Surface  Theory,”  AIAA  Journal,  Vol.  3, 
No.  1,  pp.  3-16  (Jan  1965). 

40.  Olson,  R.E.  and  W.F.  Brownell,  “Facilities  and  Research  Capabilities-High  Speed 
Phenomena  Division,  David  Taylor  Model  Basin,  Langley  Field,  Virginia,”  David  Taylor  Model 
Basin  Report  1 809  (Apr  1 964). 


